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Abstract 



To optimally utilize all the exciting cosmological data coming in 
we need to sharpen also the theoretical tools available to cosmologists. 
One such indispensible tool to understand hot big bang cosmology is 
finite temperature field theory. We review and summarise the efforts 
made by us to use finite temperature field theory to address issues of 
current interest to cosmologists. An introduction to both the real time 
and the imaginary time formalisms is provided. 

The imaginary time formalism is illustrated by applying it to under- 
stand the interesting possibilty of Late Time Phase Transitions. Recent 
observations of the space distribution of quasars indicate a very notable 
peak in space density at a redshift of 2 to 3. It is pointed out that this 
may be the result of a phase transition which has a critical temperature 
of roughly a few meV (in the cosmological units, h = c = k = 1), which 
is natural in the context of massive neutrinos. In fact, the neutrino 
masses required for quasar production and those required to solve the 
solar neutrino problem by the MSW mechanism are consistent with 
each other. As a bonus, the cosmological constant implied by this 
model may also help resolve the discrepancy between the recently mea- 
sured value of the Hubble Constant and the age of the universe. 

We illustrate the real time formalism by studying one of the most 
important time-dependent and non-equilibrium phenomena associated 
with phase transitions. The non-equilibrium dynamics of the first stage 
of the reheating process, that is dissipation via particle production is 
studied in scalar field theories. We show that a complete understand- 
ing of the mechanism of dissipation via particle production requires a 
non-perturbative resummation. We then study a Hartree approxima- 
tion and clearly exhibit dissipative effects related to particle produc- 
tion. The effect of dissipation by Goldstone bosons is studied non- 



1 



perturbatively in the large N limit in an 0(N) theory. 

We also place our work in perspective and point out some of the 
related issues which clearly need further exploration. 
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Chapter 1 



INTRODUCTION TO FINITE 
TEMPERATURE FIELD 
THEORY AND ITS 
RELEVANCE TO COSMOLOGY 



1.1 Introduction 

Cosmology is one of the most rapidly growing fields of physics today. This is 
so because there is a fairly large number of observations which we are now in 
position to make which have been impossible to do so far in mankind's history. 
The combination of both space based observations and improvements in technology 
and techniques of ground based observations is dramatically improving both the 
quantity and quality of data now available. These data will clearly lead to an 
improved understanding of the physics of the universe we live in. 

There are many examples of areas of cosmology where better data continues to 
lead to better understanding. In particular, structure fomation continues to benefit 
from data from a variety of sources. These data includes improving determination 
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of the clustering of galaxies and other larger scale structures[0]. Further light is 
being shed on structure formation by examining number densities of quasars]^], 
clustering properties of quasars and the study of Lyman-a clouds [Q. A par- 
ticularly important piece of information that has direct implications for structure 
formation has become available in recent years - the detection of anistropies in the 
Cosmic Microwave Backround Radiation(CMBR)|jl2|. This has been measured by 
COBE and continues to be pursued by ground based and balloon based experi- 
ments. We are starting to get information on the CMBR anisotropy on a variety 
of angular scales and this together with available data on the clustering of galax- 
ies and other larger scale structures is starting to severely constrain models of 
structure formation. 

Another example of fresh observations which have the potential of leading to 
an improved understanding of the physics of the universe we live in is the accu- 
rate determination of the Hubble Constant using Cepheid variables as standard 
candles[^, 0. This is among the many achievements of the new refurbished Hub- 
ble Space Telescope. The ability of the Space Telescope to see at a level of detail 
which has not been possible earlier continues to have an impact on many fields 
of astrophysics. Another of its significant contribution has been to provide com- 
pelling evidence for the existence of supermassive blackholes at the center of the 
galaxies. An examination of the rotation velocities of objects close to the center 
of the galaxy implied the eixstence of an enormously massive object confined to 
such a small region of space that there is no known mechanism of preventing the 
collapse of such a massive object to form a black hole. 

The determination of the Hubble Constant is an example where both previ- 
ously impossible space based observations^ and improvements in ground based 
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techniques and technology!^ are happening at the same time to complement and 
confirm improved, more accurate results. 

The accurate determination of the Hubble constant is also an example of a 
recent observation that has important implications for cosmology and the potential 
of leading to an improved understanding of the physics of the universe we live in. 
One can use the Hubble constant to determine the "expansion age" of the universe 
which is the age of the universe derived on the assumption that we live in a fiat 
universe with no vacuum energy. This expansion age may turn out to be less than 
the observed ages of globular star clusters. If all the data stand up to the test of 
time then clearly we have a very important piece of fresh information about the 
universe. 

One of the ways of resolving this apparent confiict between the age of the 
universe and the Hubble constant is to introduce a cosmologiacal constant in Ein- 
stein's equations which determine the evolution of the universe. The naturalness 
or unnaturalness of this solution is discussed in a later chapter and will be com- 
mented on later, but there is a further point about the role of observations in this 
context that we would like to make here. 

The existence and magnitude of the cosmological constant is an issue which 
is also benefitting from continued improvements in astronomical observations. In 
particular, since the cosmological constant modifies the geometry of the universe 
and so infiuences the propogation of light it can be observationally determined or 
at least constrained. Thus by observing such quantities as quasar lensing frequen- 
cies and comparing them to those expected for various values of the cosmological 
constant one can determine or at least constrain the value of the cosmological 
constant]^. There is admittedly a dependence on the modelling of the time evolu- 
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tion due to ongoing structure formation involved, but for a given detailed model of 
structure formation one can use the observations to shed light on the cosmological 
constant issue. Indeed, considerable progress has been made in this field too in the 
last few years and continued improvements will only make this a further source of 
rich information on cosmology. 

So far we have given examples of ongoing observations that are yielding impor- 
tant information about cosmology. There is also a host of observations that will 
be coming on-line in the future which will undoubtedly give us valuable new infor- 
mation about cosmology. In particular, neutrino astronomy]^ and gravitational 
wave astronomy^ will allow us to take a look at the universe in a way we have 
never been able to explore previously. 

It is clear that cosmology is becoming an intensely data-driven subject. The 
fresh new data that continues to stream in will undoubtably lead to an increased 
understanding of the universe we live in. However, to use all the data efficiently 
and effectively, utilize it to yield insights into the fundamental physics of the pro- 
cesses involved requires us also to equip ourselves with some theoretical tools and 
techniques. 

Field theory is clearly one of the essential tools for an understanding of cos- 
mology. In particular, finite temperature field theory plays a central role in under- 
standing hot big bang cosmology. This is perhaps underlined by the recent COBE 
observation of the perfect thermal Planckian spectrum for the photon gas that 
inhabits the universe. The Hubble expansion then tells us that at earlier times the 
energy density was higher leading to higher temperatures for the photon gas at ear- 
lier epochs. Also at the earlier higher densities the higher interaction rates would 
imply that there was thermal equilibration of the various particle distributions 
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with the photon gas at earher epochs. 

In addition to this we beheve there were a number of phase transitions during 
the evolution of the universe. These phase transitions have important imphcations 
for cosmology. One of these phase transitions is believed to be responsible for 
inflation which solves many of the problems of standard cosmology, such as the 
flatness problem and the entropy problem[||, |^. The natural tool to study all 
these phenomena is finite temperature field theory. This point has been made in 
a number of excellent contemporary books on cosmology 

In this thesis, we attempt to provide an introduction to finite temperature 
field theory geared towards those interested in cosmology. Towards this goal we 
have used the techniques of finite temperature field theory to address some of the 
subjects and issues of current interest to cosmologists. 

Phase transitions play an important role in modern theories of particle physics 
and cosmology Symmetry breaking phase transitions are an integral part of 



our current understanding of the unification of fundamental interactions 

In the Big Bang cosmology as the universe evolves it expands and it cools. 
Starting at high temperatures with a very symmetric universe as the universe 
expands and cools it goes through a sequence of phase transitions. Thus, first the 
GUT symmetry occurs followed by the electroweak sjTiimetry breaking followed 
by other phase transitions. 

Among the interesting cosmological issues that are intimately linked with these 
phase transitions are the formation and evolution of topological defects PUl E3 



|3, ITB[, infiation[iO, 53, H^, EO], structure formation and baryogenesis. 



It is clear that we need a systematic framework in which to address all of 
these issues. Finite temperature field theory pD|, ^ Q ^ is definitely the 
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framework one needs to use to get detailed quantitative information about these 
interesting and important events. Within the framework of finite temperature field 
theory there are two different formalisms one could use. These are respectively the 
imaginary time formalism and the real time formalism of finite temperature field 
theory. Depending on the particular questions or issues one wishes to address one 
formalism may be better than the other and on occasion one of the formalisms 
may not be relevant to the issues to be addressed. In particular, to determine the 
existence of phase transitions, the magnitude of the critical temperature and the 
nature of topological defects formed at the phase transition one may wish to use 
the imaginary time formalism. However to study the time dependent and non- 
equilibrium phenomena associated with phase transitions the natural tool to use 
is the real time formalism of finite temperature field theory. The various uses of 
the two different formalisms of finite temperature field theory will be illustrated in 
detail by studying specific examples of current interest. 

Phase transitions that occur after the decoupling of matter and radiation have 
been discussed in the literature as Late Time Phase Transitions (LTPT's). The 
original motivation for considering LTPT's|ll| [|l^ |l7| |25| [20| p3|was the need 



to reconcile the extreme isotropy of the Cosmic Microwave Background Radiation 



(CMBR)[0 with the existence of large scale structure 1 13] and also the existence 
of quasars at high redshifts |I^] . 

Discussions of realistic particle physics models capable of generating LTPT's 



have been carried out by several authors 1 20 1 pSl. It has been pointed out that the 



most natural class of models in which to realise the idea of LTPT's are models 
of neutrino masses with Pseudo Nambu Goldstone Bosons (PNGB's). The reason 
for this is that the mass scales associated with such models can be related to the 
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neutrino masses, while any tuning that needs to be done is protected from radiative 



corrections by the symmetry that gave rise to the Nambu-Goldstone modes|21 



Holman and Singh p3[] studied the finite temperature behaviour of the see-saw 
model of neutrino masses and found phase transitions in this model which result 
in the formation of topological defects. In fact, the critical temperature in this 
model is naturally linked to the neutrino masses. 

The original motivation for studying the finite temperature behaviour of the 
see-saw model of neutrino masses came from a desire to find realistic particle 
physics models for Late Time Phase Transitions. It now appears that this may 
also provide a physically appealing and observationally desirable magnitude for 
the cosmological constant. 

In particle physics one of the standard ways of generating neutrino masses has 



been the see-saw mechanism ||22|. These models involve leptons and Higgs fields 
interacting by a Yukawa type interaction. We computed the finite temperature 
effective potential of the Higgs fields in this model. An examination of the manifold 
of degenerate vacua at different temperatures allowed us to describe the phase 
transition and the nature of the topological defects formed. 

To investigate in detail the finite temperature behaviour of the see-saw model 
we selected a very specific and extremely simplified version of the general see- 
saw model. However, we expect some of the qualitative features displayed by our 
specific simplified model to be at least as rich as those present in more complicated 
models. 

In the following chapter we will study in some detail a specific simple version 



of the see-saw model of neutrino masses|22]. By studying the finite temperature 



behavior of this model Holman and Singh demonstrated the existence of phase 
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transitions, establishing the magnitude of the critical temperature. This was done 
by using the imaginary time formalism of finite temperature field theory. By study- 
ing the manifold of degenerate vacua we also established the nature of topological 
defects that would be formed. The phase transition implied by this model may 
have interesting cosmological consequences if we use the neutrino masses implied 
by the MSW solution||2^ to the solar neutrino problem. In particular, these phase 
transitions might help explain the peak in the quasar space densities observed at 



redshifts between 2 and 3[|T6l. Also, the magnitude of the cosmological constant 
implied by this model might help resolve the apparent conflict between the age of 
the universe and the Hubble constant |^]. We will study the details of the finite 



temperature see-saw model in Chapter 2. 



The real time formalism of finite temperature field theory pSj ^ ^ has been 



used by us to study time dependent and non-equilibrium phenomena during and 
following a phase transition such as domain formation and growth |T^, thermal 



activation ||60||, dissipation via particle production and reheating ||6T1, |62| . In Chapter 



3 we will study in detail the mechanism of dissipation via particle production as 
an illustration of the use of the real time formalism to study time dependent 
phenomena following a phase transition. 

We next review the formalism we will use in later chapters. The remainder of 
this chapter will be devoted to introducing first the imaginary and then the real 
time formalisms of finite temperature field theory. 



14 



1.2 The Imaginary Time Formalism 



In finite temperature field theory the basic quantity of interest are the finite tem- 
perature Green's functions. To study phase transitions and identify the critical 
temperature and nature of topological defects one uses the effective potential which 
can be expressed in terms of 1 particle irreducible Greens functions as discussed 
below. 

Let us quickly review the zero temperature field theory and the zero tempera- 
ture effective potential first and then discuss their finite temperature counterparts. 
We will consider the simplest possible example of a scalar field theory to introduce 



the formalism 39 



The starting point of our discussion is the generating functional for the Greens 
functions, 



W[J] = J D(j)exp |z J d^x{L + J0)| 



where L is the Lagrangian density describing the field theory of interest, (p is the 
scalar field of interest and J is the source which allows us to express the n-point 
greens functions in terms of the generating functional in the following way: 



,X2,...,Xn) - -ny^r 



6J{Xi)...6J{Xn) 



1.21 



J=0 



Similarly we can introduce the quantity, Z[J] = —i lnW[J], which can be 



shown to be the generating functional of connected greens functions p9 



X2, 



(5"ln^[J] 

6J{Xi)...SJ{Xr, 



1.3) 



j=o 
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Then, the classical field which is the vacuum expectation value of the field is 
given by, 



(o|o) 5j{xy ^ > 

The effective action, F is the functional Legendre transform of the generating 
functional Z[J], 

r[0,] = Z[J\ - I d^xJ{x)(t>,{x). (1.5) 
One can also expand T[(j)^ in powers of 0c 

r[0c] = E~T / d:^Xi...d}XnT^"\xi,X2,...Xn)(pc{Xi)...(pc{Xn)- (1-6) 

It can be shown that, r(")(xi, X2, ...x„) = (xi)''^54> (x ) -'■ P^-^ticle ir- 

reducible (IPI) n-point greens functions. Recall that the IPI Greens functions 
receive contributions only from Feynman diagrams which cannot be disconnected 
by cutting any one internal line. 

Making a derivative expansion of the effective action allows us to isolate the 
effective potential as the lowest order term in the expansion. 



r[0e] = / d'x [-K=//(0c) + \{d^ct>cfF{<t>,) + 



;i.7) 



In particular, for (pci^x) = (f>c = constant, we have (9^0c = 0. Thus, by taking 
Fourier transforms one obtains from the previous equations the following relation. 
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1.2.1 Calculation of the effective potential using the tad- 
pole method 

The effective potential formalism is the most well-known and widely used formal- 
ism to study the phenomenon of spontaneous symmetry breaking in field theory. 
The formalism is by now well-developed and there exist a number of excellent dis- 



cussions on the effective potential and methods to calculate it[|J, The most 
straightforward and frequently discussed method of calculating the effective po- 
tential is by evaluating all the IPI n-point Greens function and then explicitly 
performing the sum over n as required by the relation given earlier to obtain an 
expression for the effective potential. There exists another method which is less 
well-known, a little devious, usually easier to implement in practice and is our 
method of choice to calculate the effective potential. By being a little devious 
we can express the effective potential completly in terms of a 1-point IPI Greens 
functions and this makes the computaion of effective potential much simpler in 
practice. We now turn to a description of this method and its illustration using a 
simple example. 

By introducing a shift in the fields, 

(1.9) 

one can obtain an expression for the effective potential in terms of the IPI n-point 
functions in the shifted theory, F^") which is more useful for actually calculating 
the effective potential in practice. Thus, 



K//(0c) = - E = 0, = 0) (0, - ujY. (1.10) 



This implies. 



17 



= -rL^Hpi = o). (1.11) 

Therefore, we have, 

K// = -/r«da;. (1.12) 

Note that by making the above transformations we can calculate the effective 
potential by knowing only the 1-point 1 PI function in the shifted theory as opposed 
to calculating all the n-point functions and then summing over n in the original 
expression for the effective potential. This method is called the tadpole method 
because the calculation of the IPX 1-point function, F^^^ involves the computation 
of a Feynman diagram which bears a closer resemblance to a tadpole than the 
penguin diagram's resemblance to a penguin. 

Let us illustrate the method with a simple example. Consider a scalar field 
theory with the following tree level potential, 

= Im'cP' + (1.13) 
Performing the shift, (p ^ (p — ui , we get, 

^ lm\<l> - + ^{<l> - u;)\ (1.14) 

In the shifted theory the mass of the field, determined from the coefficient 
of the 0^ term is given by, M^{u) = w? + The computation of F£^' involves 
the cubic self-coupling i.e. the coefficient of the 0^ term which is ^uj. Thus, the 1 
loop contribution to F£^^ is, 
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d(f) 



dVeff 



(1.15) 



dcu 



Thus, 



eff 



In 



k^-Mj{uj)+ie 

A2 



:i.i6) 



By performing the rotation to Euchdean space the above integral can be eval- 
uated to yield the 1 loop contribution to the effective potential for this theory to 
be, 



A' 



eff 



647r2 



In 



A2 



(1.17) 



where A is the momentum cutoff introduced to regularise the integral. This theory 
can now be renormalized by adding counterterms and specifying a renormalization 
prescription as is done below. The counterterms we add are 



;i.i8) 



The total potential, V then becomes the sum of the bare potential and SVct-We 
specify the following renormalization prescription: 



J 4>=o 



d^V 
1^ 



Xr , y(0 = 0) = Vo. 



;i.i9) 



J </.=0 
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This yields the following expressions for the renormalized parameters: 




Vo 



Ay 



A + a 



5m^ + 



6V + 



327r2 
AA 



4 

m 



327r2 



327r2 

3A2 r 




(1.22) 



(1.20) 



(1.21) 



(1.23) 



Thus in terms of these renormalized parameters the complete effective potential 
becomes, 



where m| = + 

In equilibrium finite temperature field theory we are interested in evaluating 
the finite temperature Green's functions, 



Recall that P — where T is the temperature. The temperature is a well 
defined quantity only for systems in equilibrium. In cosmology, we are frequently 
interested in systems which are in local thermal equilibrium. This happens when 
the interaction rates are high compared to the expansion rate of the universe. In 
such situations the above quantity is meaningful and of interest. There are of 
course, other situations such as when a phase transition occurs that the system is 
far from equilibrium and we will later discuss a formalism which allows us to keep 
track of the time evolution under these circumstance. For now, let us concentrate 





Gt (-^1) -^2, ■■■Xn) 




(1.25) 
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on equilibrium finite temperature field theory. The fundamental quantity of inter- 
est for equilibrium finite temperature field theory is the partition function which 
is the generating function of Green's functions, 



(1.26) 



(1.27) 



Thus we need to evaluate (0|e ^^|</>). Recall that similar quantities, 



(g';t'|g;t) = (g'|e-^(*-*')|g) 



(1.28) 



are frequently evaluated in quantum mechanics. This points the way to using 
well developed techniques in zero temperature field theory for calculations in finite 
temperature field theory by making the identification ir = — /3. Thus the inverse 
temperature is formally identified with the imaginary time direction and it is for 
this reason that the formalism described in this section is called the imaginary 
time formalism of finite temperature field theory. 

An important point to notice is that since we are interested in evaluating the 
trace, the ' initial ' and ' final ' states have to be identified leading to periodicity 
in imaginary time r for bosons. Fermions are anti-symmetric under the exchange 
of states and this leads to an anti-periodicity in imaginary time for the fermions. 

Knowing the Hamiltonian allows one to extract the finite temperature Feynman 
rules which are then used to calculate the finite temperature effective potential 
using the tadpole method. Let us do this in some detail for the X(f)'^ theory. The 
most important modification to the Feynman rules at finite temperature is in the 
propogator. In fact, the factors at the verices and the combinatorics are the same 
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at finite temperature as they are at zero temperature. The important modification 
of the propogator comes essentially from the fact that we have definite periodicity 
in imaginary time for the bosons and anti-periodicity in imaginary time for the 
fermions. To determine the propogator it is sufficient to consider only the quadratic 
part of the action, 



Jo 



drd X c 



Introducing the Fourier expansion for the field, 



0. (1.29) 



-1 oo „ 

H^,^) = g E j 4>m{k)e'^-^+^'"'^^d'k (1.30) 



/3 

one gets the propogator in momentum space to be. 



1 



D{k,u:,m) = ^^-TT-TT. (1-31) 



where the quantity ui2m = results from the periodicity in imaginary time. Also 
notice that the periodicity in imaginary time leads to a dicrete sum in the fourier 
transform. Thus, 

For fermions we have anti-periodicity in imaginary time and thus instead of 
uj2mi the quantity entering expressions is, uj2m+i — ^^^2l!^^. 

Once again we will use the tadpole method, this time to compute the finite 
temperature correcions to the effective potential. The relation which gives the 
finite temperature correction, AVr is, 
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dAVr 1^ f d^k Xoo 



(1.34) 



duj P^J {2nful^ + k^ + Ml{uj) 



I f d^k ^ 
= 7522 77^^^ 



A2 



So far, the only thing that is really different from the zero temperature discus- 
sion has to do with the fact that there is a periodicity in imaginary time leading to 
a discrete fourier sum in the timelike direction. To proceed further, we choose to 
evaluate the sum first and then evaluate the integral. The sum can be evaluated 
by using the identity, 

oo 

ln(a;2„ + E^) = (3Ek + 2 ln(l - e"^-^*) + {Ek independent terms) (1.35) 

m=— oo 

where 

El^e + Ml,u;,m-^. (1.36) 

Using the above identity we can evaluate the finite temperature correction to 
the effective potential. 



f d^k , 


1 - e-^^^" 


/ (2.)3 





(1.37) 



At this point we note that the second term on the right is finite and the first 
term on the right hand side contians all the divergent part. Further, wc note that 
in fact the divergence structure of this first term is identical to the divergence 
structure of the zero-temperature part of the effective potential. It thus follows 
that the counterterms that we introduced to absorb the divergences of the zero- 
temperature effective potential will be sufficient to absorb all the divergences of 
the finite temperature corrections to the effective potential. Hence, the net effect 
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of the finite temperature corrections is to add only tfie second term on the right 
hand side (which we wiU caU /S.Vj< to the expression for the renormahsed effective 
potential. Let us evaluate the consequences of this term in some further detail. 
First, let us set x — 13k and do the solid angle integral. This yields, 

1 



dx x^ In 



-(a;2+/32M2(</,e))2 



(1.38) 



For the region where << 1 we can do a high temperature expansion to the 
above integral which yields. 



90 



647r 



24 



1271 



2^2 



+ 



(3/2-21n47r-27) 



;i.39) 



To lowest order. 



(1.40) 



90 24 

Thus, if we collect the coefficient of the 0^ from the zero temperature potential 
and the finite temperature correction and define this coefficient as w?{T), then we 
get. 



m^{T) = m|=o + 
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;i.41) 



which indicates how a theory which has a broken symmetry at zero temperature 
(m|,^Q < 0) can have the symmetry restored at some higher temperature such 



that w?{T) — m\^Q + ^ > 0. Thus, the critical temperature will be given by 
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m\T,) = = ml^^ + ^. 

Thus we see how the imaginary time formalism of finite temperature field theory 
can be used to determine the critical temperature given a model of the interactions 
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of fields we are interested in, in terms of the parameters of the model. In this section 
we have used the simplest possible model of interactions to introduce the concepts, 
formalism and techniques of the imaginary time formalism of finite temperature 
field theory. In chapter 2 we will use this same formalism and these same techniques 
to study the finite temperature behaviour of a specific see-saw model of neutrino 
masses. As we will see the finite temperature behaviour of this model may have 
many interesting cosmological consequences including a possible explanation of 
the peak in comoving quasar space densities as well as providing a cosmological 
constant of the right magnitude to remove the discrepancy between the Hubble 
constant and the age of the universe. Of course, in addition to being of current 
interest this model will also allow us to illustrate the techniques and utility of the 
imaginary time formalism of finite temperature field theory. In particular, we will 
be able to analyze the manifold of degenerate vacua of the effective potential in 
that model to extract information about the nature of topological defects that will 
be formed at various temperatures. 

1.3 The Real Time Formalism 

So far we have studied the imaginary time formalism of finite temperature field 
theory. As already pointed out however, the imaginary time formalism cannot be 
used to study time dependent phenomena since the ' imaginary time ' is identified 
with temperature. The best one can hope to do in this formalism is to stipulate 
that the temperature is some specified function of time T{t). However, this implies 
that the system is always very close to some equilibrium configuration with a well 
defined temperature T. In particular, this implies that the density matrix is always 
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expressible in the form p — e~^^ . As we shall see during the transition period 
immediately following a phase transition the distribution function does not look 
like it can be fit by a purely thermal distribution function with a single well defined 
temperature at each instant of time. 

The physical quantity best suited to describe the time evolution of a quantum 
statistical system out of equilibrium is the density matrix, p{t). The density ma- 
trix can be expressed in terms of a complete set of eigenstates in the following 
way, p = J2i \'^i)Pi{'^i\ where pi represents the probability of finding the state of 
the system in the state Note that the density matrix need not neccesarily 

be diagonal, however it must be hermitian and hence can be diagonalised by a 
unitary transformation and expressed in the form given above. The density ma- 
trix is particularly useful for describing the state of a system which is not a pure 
quantum mechanical state but is instead in a mixed state. In such cases the only 
information we may possess is the statistical information contained in the probabil- 
ities of finding the various different eigenstates of the system. An example of such 
a mixed statistical system is completely unpolarised light for which the density 
matrix would be an equally weighted sum over the different polarization states of 
light. Another example of a mixed statistical system which occurs frequently in 
nature and will be of concern to us is a system in thermal equilibrium at a given 

temperature T. For such a system all that we may know is that the probabihty of 

-El 

finding an eigenstate with energy Ej is proportional to e ^t. For statistical sys- 
tems we are frequently interested in the time evolution of a specific initial mixed 
state under the influence of a time dependent Hamiltonian. 

The time evolution of the system is then determined by the quantum Liouville 
equation, 
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ih^^ [Hit), Pit)] (1.42) 

where Hit) is the Hamiltonian of the system. 

The formal solution of the quantum Lioville equation is given by, 

pit)^Uit,to)pito)U-\t,to) (1.43) 

where p(io) is the intial density matrix at time to and U it, to) is the time evolution 
operator from to to t. 

Typically wc are interested in systems that are initially in thermodynamic 
equilibrium with a well defined initial inverse temperature /?j = "^^^ initial 

density matrix is then given by, 

p(^o)=e-^'^^ (1.44) 

Notice that if wc arc interested in the ground state of the system which is the 
equilibrium state at zero temperature, we can isolate it by taking the limit /?j — > oo. 

The initial Hamiltonian Hi which governs the behavior of the system for times 
t ^ to prepares the system to be in equilibrium with a well defined temperature and 
Pi uptil to- Interesting non-equilibrium phenomena will occur if for times t > to, 
the Hamiltonian becomes H^^oiit) which is different from the initial Hamiltonian 
Hi. The time evolution of the expectaion value of any operator O is given by, 

Tr[e-P^^^U-\t)OUit)] 
Tr[e-^.H^] ■ 

The above expression can be put into a more compact form by noting once 
again that the time evolution operator, U also involves the exponential of the 
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Hamiltonian. In particular, for any time T < to, 



U{T) = exp[~iTHi]. (1.46) 

We can use this fact to absorb the exp[—f3iHi] factor into a modified time 
evolution operator. Let us re-express, exp[— = exp[—iHi(T — i[3i — T)] = 
U{T — i(3i,T). Further, we insert the factor U^^{T)U{T) = 1 in the numerator 
and using the fact that U^^{T) commutes with exp[— we can re-express, 

_ Tr[UiT-^p.,t)OU{t,T)] 

Tr[U{T-z(3,,t)] ■ ^'-^'^ 

We have used slightly convoluted means to arrive at a compact result. In 
particular, we have hidden the initial density matrix but we have entered the 
complex time plane. Notice that in the imaginary time formalism we had only 
the imaginary time direction playing any role. In the real time formalism both the 
imaginary direction as well as the real direction enter the discussion. Thus the real 
time formalism should really be called the complex time formalism. However, that 
might have made the real time formalism sound even more dauntingly complex 
than it already appears. This may have impeded progress in this field for a few 
more years in real time. 

We will use the real time formalism to study in detail the mechanism of dissipa- 



tion via particle production [^T], ^ in a later chapter. This is a process that is very 
important for understanding reheating following inflation. In chapter 3 we use the 
real time formalism to study in detail one of the most important time dependent 
and non-equilibrium phenomena associated with phase transitions. This is the 
phenomenon of dissipation of the initial vacuum energy of the field. The physical 
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mechanism for this is the particle production which takes place as the field oscil- 
lates about the minimum of its potential. Wc used the quantum Lioville equation 
to study the time evolution of the system starting with an initial configuration 
described by a density matrix, Piniuai- 

Once again the motivation for presenting this discussion in detail is twofold. 
First, it is a subject of current cosmological interest in its own right. Secondly it 
provides a detailed illustration of the techniques that will undoubtably be useful in 
the study of a number of other time-dependent and non-equilibrium phenomena. 
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Chapter 2 

IMAGINARY TIME 
APPLICATION: A LATE TIME 
PHASE TRANSITION LINKED 
WITH MASSIVE NEUTRINOS 



2.1 Introduction 

One application of the imaginary time formalism that is of current cosmological 
interest is to the treatment of Late Time Phase Transitions(LTPT's). We will first 
describe what LTPT's are and then point out why we think they have a number of 
extremely interesting and potentially important cosmological consequences. Once 
this is done we will turn to a detailed discussion of a compelling particle physics 
model for such phase transitions. The detailed investigation of this model will also 
illustrate some of the strenghts of the imaginary time formalism for determining 
the existence of phase transitions, the critical temperature of the phase transition, 
the existence and nature of topological defects formed during the phase transition 
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in a given model for interacting fields. 

Phase transitions that occur after the decoupling of matter and radiation have 
been discussed in the literature as Late Time Phase Transitions. The original 



motivation for considering LTPT's|TT|, [T6|, |T^, |T9|, ^ p^was the need to reconcile 



the extreme isotropy of the Cosmic Microwave Background Radiation (CMBR)p!2| 
with the existence of large scale structure and also the existence of quasars at 
high redshiftsin]. The hope that LTPT's would not disturb the CMBR was not 



realised[15|. In light of the fact that the anisotropy measurements of the CMBR 



are getting better the distortions of the CMBR predicted by LTPT's may 
become important observational tests of these models. 

In our mind the current motivation for a continued interest in LTPT's is 
twofold. First, LTPT's may provide an explanation for the peak observed in co- 
moving quasar space densities plotted as a function of redshift . As we will see 
in detail in the discussion later, the most compelling particle physics model for 
LTPT's is the finite temperature behaviour of a see-saw model of neutrino masses. 
In such a model one finds a critical temperature which is a few times the relevant 
light neutrino mass. If we input into this model the neutrino masses required to 
solve the solar neutrino model by the MSW mechanism Q we obtain that the 
phase transition happens at the right epoch to explain the peak in quasar space 
densities. A second important motivation for a continued interest in the LTPT 
model being presented in this chapter is that it may help remove the discrepancy 
between the recently determination of the Hubble constant and the age of the 
universe. 

LTPT's are used to generate density fluctuations either by bubble nucleation 



dynamics of a slow rolling field or by the formation of topological defects 
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such as domain walls(HSF). The appeal of LTPT's for structure formation hes in 
the fact that the density perturbations created are immediately non-linear and so 
will lead to much faster formation of structure. The models we will be considering 
in detail in this chapter involve the formation of topological defects as a result of 
LTPT's. The original idea of the production of soft domain walls as a result of an 
LTPT (HSF) has received several criticisms. A problem with the use of domain 



walls to form Large Scale Structure has been pointed out by Kolb and Wang[|I| 
They argue that domain walls may never form in a LTPT since thermal effects 
may not be able to drive different regions of the Universe to different parts of the 
disconnected vacuum manifold. However they also point out that domain walls 
would form if large fluctuations in the scalar field exist before the phase transition. 
One of the main critcisms of soft domain walls is that they may produce significant 
distortions of the CMBR. The most significant microwave distortion comes from 
collapsing domain wall bubbles [P!^] . These distortions produce hot and cold spots 
on ~ 1° angular scales and provide a signature for the formation of domain walls 
after recombination. 

Since the collapse of closed domain wall bubbles to form black holes which 
then act as the central engine of quasars is one of the main themes of this part of 
this thesis, we now turn to a quick discussion of the distortions of the CMBR this 
produces. The temperature shift due to a photon traversing a collapsing domain 
wall bubble is 

AT 



T 



2.64 X 10-%'^(3Aa/{10MeV^) (2.1) 



where h, A, j3 are dimensionless numerical constants of order unity and a is the 
surface tension of the domain wall. The present measurements of the CMBR 
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anisotropy then imply|^^ that a < 0.5 MeV^. This constraint, though important 
to keep in mind is not a problem for the viabihty of the model being presented 
in this chapter. We will return to this point after discussing our particle physics 
model for the LTPT. 



2.2 Particle physics models for LTPT's 

Discussions of realistic particle physics models capable of generating LTPT's have 
been carried out by several authors[|20|, It has been pointed out that the most 
natural class of models in which to realise the idea of LTPT's are models of neutrino 
masses with Pseudo Nambu Goldstone Bosons (PNGB's). The reason for this is 
that the mass scales associated with such models can be related to the neutrino 
masses, while any tuning that needs to be done can be protected from radiative 



corrections by the symmetry that gave rise to the Nambu-Goldstone modes|21 



Holman and Singh p3[] studied the finite temperature behaviour of the see-saw 



model of neutrino masses and found phase transitions in this model which 
result in the formation of topological defects. In fact, the critical temperature in 
this model is naturally linked to the neutrino masses. 

To investigate in detail the finite temperature behaviour of models allowing 
for natural LTPT's we selected a very specific and extremely simplified version of 
the general see-saw model. However, we expect some of the qualitative features 
displayed by our specific simplified model to be at least as rich as those present in 
more complicated models. 

We chose to study the 2-family neutrino model. Because of the mass hierarchy 
and small neutrino mixings]^] the expectation was to capture some of the essential 
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physics of the z/g-z/^ system in this way. The 2-family see-saw model we consider re- 
quires 2 right handed neutrinos N}^ which transform as the fundamental of a global 
SUr{2) symmetry. This symmetry is implemented in the right handed Majorana 
mass term by the introduction of a Higgs field aij, transforming as a symmetric 
rank 2 tensor under SUpi{2) (both and cjij are singlets under the standard 
model gauge group). The spontaneous breaking of SUr{2) via the vacuum ex- 
pectation value (VEV) of a gives rise to the large right handed Majorana masses 
required for the see-saw mechanism to work. Also, the spontaneous breaking of 
SUr{2) to U{1) gives rise to 2 Nambu Goldstone Bosons. The SUr{2) symme- 
try is explicitly broken in the Dirac sector of the neutrino mass matrix, since the 
standard lepton doublets II and the Higgs doublet $ are singlets under SUr{2). It 
is this explicit breaking that gives rise to the potential for the Nambu Goldstone 
modes via radiative corrections due to fermion loops. Thus, these modes become 
Pseudo Nambu Goldstone Bosons (PNGB's). 

The relevant Yukawa couplings in the leptonic sector are: 

- C-yu^ = yaMNii^ + yNj^Ni "aij + h.c. (2.2) 

where a,i,j = 1,2. The S'f/ij(2) symmetry is implemented as follows: 

cTi, ^ U^MU^Yj (2.3) 

where f/j is an SUfi{2) matrix. The first (Dirac) term breaks the symmetry ex- 
plicitly. 

We can choose the VEV of cr to take the form|^: (cTjj) = f6ij, thus break- 
ing SUii{2) spontaneously down to the f/(l) generated by T2 (where are the 
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Pauli matrices). This symmetry breaking gives rise to the 2 PNGB's, whose finite 
temperature effective potential is of interest. 

After the Higgs doublet acquires its VEV, we have the following mass terms 
for the neutrino fields: 

- JC^sss = rUaiPlNl, + MNrUU^N% + h.c. (2.4) 

where are the standard neutrinos, niai = Uai 'y/v^, M = yf/\/2. 

Examining the expression for — vCmass above we notice that it may be expressed 
more compactly by introducing the fields, 

N'r = U'^Nr (2.5) 
Then the expression for — £mass becomes, 

- >C^ass = rUaiPlUiN'i + MN'j^N'i + h.c. (2.6) 

At this point we can introduce the u - mass matrix to rewrite the above terms 
in a block matrix form : 



where the mass matrix, Alj,, 



~ I {mUf M J 
can be re-expressed in the suggestive form. 



) (2.8) 



= n ^^] + { ^^rnr n (2-9) 



Mj ' \(mUy 
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In particular, in the physical case where |maj| << M we can use perturbation 
theory techniques to diagonalisc the mass matrix. Specifically, we need to find the 
eigenvalues of the mass matrix by using the second matrix as a perturbation on 
the first matrix in the expression above. Thus, the zeroth order eigenvalues are 
0, 0, M, M and the corresponding eigenvectors can be taken to be of the block form 
( ( ^\ ^ ^ ^ where 

e.= (;),e.= («) (2.10) 



oy Voy' Vei/' Ve2 



The most significant effects due to higher order corrections will clearly be to 
the first two eigenvalues which are both at the zeroth order in perturbation 
theory. Let us calculate the corrections to these two eigenvalues explicitly. Call 
the perturbation matrix Mi i. e. 



Then, the first order corrections to the eigenvalues A, for i = 1, 2 is given by. 



Af) = {i\Mi\i) (2.12) 

= T)(o) 

= (2.14) 

Thus the first non-zero correction to these eigenvalues occur at second order in 
perturbation theory which is given by. 



- 1^ 7% y^-^^) 

3+i \ ~ ^3 
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Substituting the explicit form for Mi and the eigenvectors in the above expres- 
sion we get, 



1 



{mU)l^ + (mU) 



12 



(2.16) 



and 



A? 



M 



{mU)l^ + (mU) 



22 



(2.17) 



where 



(mU)aj = rUaiUi (2.18) 
= ma^{exp{^^■T/fy, (2.19) 

For the case we are interested here U belongs to SUr{2) and can be expressed 
in the convenient form, 



U = exp[ze-T//] (2.20) 
= cos||e||+if-Tsin||e|| (2.21) 

where | |e| | = + ^| + and = 6/ V^f+^l+l- 

Recall that wc chose the symmetry breaking scheme where SUr{2) was spon- 
taneously broken down to the U{1) generated by T2. Thus the Nambu-Goldstone 
modes point in the t\ and T3 directions. The most general form of U that we will 
therefore be interested in therefore is of the form. 



C/ = exp[i(eiTi+e3T3)//] 
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(2.22) 



In particular, 



cos||e|| + ^|3sm||C|| 
^fisinllell 

cos 1 1^1 1 -ifssin 



If for concreteness and simplicity we take niai — m Sai then, 



Al- 



es) 



m 



cosmiD+^sMmw) 



(2.23) 
(2.24) 
(2.25) 
(2.26) 



(2.27) 



and 



,(2) _ ..^ 



m 



cos(2||e||)-i|3sin(2||e||) 



(2.28) 



Since the effective potential for the ^j's is generated by the propogation of neu- 
trinos in the loops the effect of the heavier neutrinos is supressed by extra powers 
olm/M. Thus, it is the two lightest neutrinos that contribute most significantly 
to the effective potential for the .^'s. 

Substituting the above eigenvalues back in the expression for — >Cyuk and recall- 
ing that. 



;i+75)^ 



(2.29) 
(2.30) 



we get that. 
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2 

->Cyuk = -— ^i[cos(2| 1^11)+ i6sm(2| 1^1 1)75] V'l 

2 

777^ r ^ -1 

- M ^2 [cos(2| I) - sin(2| 1)75] (2.31) 

2.3 The effective potential for our model and its 
implications 

We will treat the ^^'s as classical background fields and not allow them to propogate 
in loops. Thus, the effective potential of the ^j's is going to be due to the fermion 
loops. Thus when we perform the shift ^ ^ ^ — w we need only keep the lowest 
order term in ^. (The terms that will contribute to the tadpole diagram have to 
be of the form ip^ip and can be isolated from the above expression for i2yuk by 
performing the shift C ^ C ~ ) 

When we examine the expression for >Cyuk above we notice that on performing 
the shift ^ — > ^ — w we will need to expand the terms ||^ — w|| and ]]|^E^^- To 
linear order in ^, their expansions are, 

\\^-w\\ = \\nj\\-^ + 0{e) (2.32) 

I lu"! I 

— = -w^ + —— - w^ijT—r^ + ) (2.33) 

116 -^sll \\w\\ IfII 

Using these expansions in >Cyuk and performing straightforward trigonometric 
and algebraic manipulations one can extract the following ■^^■^ couplings that are 
all relevant for computing the tadpole diagrams that will give us an expression 
for the effective potential. The coefficient of the V'l^V'i term is: — ^ [z/i + ^75^/2], 
where, 
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y2 



wl cos(2 


^11) 




\w\ 


2 



+ 



\w\ 



\w\ 



sm(2||w| I) 



(2.35) 



The coefficient of the ipi^iipi term is: — ^ [ys + nbUi], where, 



Z/3 = 



2wi sin (2 


HI) 




\w\ 





W3W1 cos(2 


\w\\) 






2 


\w\ 


3 



sin(2||w||) 



(2.36) 
(2.37) 



The coefficient of the ip2^i'4^2 term is — ^ [7/3 — ^751/4] and the coefficient of the 
'026'02 term is [yi - n5y2]- 

The tadpole diagram with one external ^1 leg will receive contributions from 
the loops due to ipi and ips and is given by, 



eff 



dwi 



Wl,W3 



= -8i 



J (2^ 



(27r)4 M F - {Al - Bl) 
where ya, were given earlier and Ai, Si are given by: 



(2.38) 
(2.39) 













Si = 


irn? 


M 





W3 



\w\ 



sin(2||w| 



(2.40) 
(2.41) 



Similarly the tadpole diagram with one external ^3 leg will be given by. 
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The above expressions can be integrated to yield, 



where, 



A4^ = ^(cos^2||e||+8sin22||e||) (2.45) 

The above expression can be exphcitly evaluated by performing a rotation to 
Euchdean space to yield. 



K//(6,6 



4 



647r2 



2 \ A? 



(2.46) 
(2.47) 



where as always A is the cut-off introduced to regularise the integral. 

In order to look for phase transitions in this system we should now add in the 
effects of the heat bath. This should be done with some caution. The reason for 
this is that by the time of decoupling, the neutrinos have long since dropped out 
of thermal equilibrium with the ambient heat bath, which would seem to indicate 
that the finite temperature formalism is inapplicable. However, we should recall 
that despite being out of thermal equilibrium, neutrinos preserve a thermal-like 
distribution after they decouple (just as the microwave background photons do). 
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with a temperature inversely proportional to the scale factor Thus, we can 
still compute the finite tempereature effective potential and use it as a diagnostic 
of critical behavior, as long as we recall that the relevant temperature is that of 
the neutrinos and not that of the heat bath. Another caveat is that once the 
temperature is comparable to the mass of the neutrino, the distribution function 
will cease to look thermal. 

The finite temperature effects can be studied by calculating the finite temper- 
ature corrections to the effective potential which is given by. 



AVt = 7-^ / dx In 

ZTT^ JO 



1 + exp[-Wx2 + -j^ 



(2.48) 



Recall that we parametrized the PNGB's as ^1 and ^3 and defined the quantity 
M to be: 

-M' = ^(cos2 2||e||+llsin2 2||e||) (2.49) 

Performing the high temperature expansion of the complete potential and discard- 
ing terms of order {M.'^Y /T"^ or higher we get. 



'^^ - + (ml + TyQ)M' + ^-^{n - log 

9U OTT^ 



Kot(6, ^3) = ViM') = (Vo - ^^) + (mt + T'/Q)M' + ^^{n - log 

(2.50) 

where n = 27 — 1 — 2 log tt ~ —2.1303, is a parameter in the model and fi is the 
renormalisation scale. A4 is naturally of the neutrino mass scale in this model. 

A study of the manifold of degenerate vacua of the effective potential at dif- 
ferent temperatures revealed phase transitions in this model accompanied by the 
formation of topological defects at a temperature of a few times the relevant neu- 
trino mass . Typically at higher temperatures the manifold of degenerate vacua 
consisted of a set of disconnected points whereas at lower temperatures the man- 
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if old was a set of connected circles. Thus, domain walls would form at higher 
temperatures which would evolve into cosmic strings at lower temperatures. 

An important point is that the potential depends on the ^'s only through A4'^. 
This has the effect of taking a two dimensional problem and turning it into a one 
dimensional one, i. e. we need only find the value of Ai^ that minimizes the 
effective potential. A picture of the potential for ml = 1 is given in fig. 1. 

We can understand some features of the manifold of degenerate vacuua of this 
effective theory if we first uncover some of the symmetries of the potential. Note 
that is left invariant by the transformations: 

iieii-iie'ii = iieii + ^ 

iieii iieii iieii ^ ■ ' 

where n is an arbitrary integer. Furthermore if we define = cos^2||^|| + 
^3 sin^ 2| |,^| I, we see that Ai"^ is restricted to lie between and 1. It vanishes when 
^3 = 0, ^1 = {2k + l)7r/4 , while it attains its maximum value of 1 when either 
^1 = or 1 1^1 1 = k7r/2 {k e Z). Thus, if = is the ground state, then the 
vacuum manifold consists of a discrete set of points, while if Ai^ = 1 minimizes 
the potential, then this manifold is the set theoretic union of an infinite number of 
circles with discrete radii and the fine = 0, which intersects all of these circles. 
We should note that due to the above symmetries some of these vacuua will become 
identified. For example, in the A4^ — case , the only distinct vacuua are: ^3 = 
and ^1 = {2k + l)7r/4, with A; = 0, 1, 2, 3. For the — 1 situation the distinct 
vacuua are = or 1 1^| | = k7i/2 , k = 0,1,2, 3. 

At this point we can define rrii, — m^/M and scale all dimensionful quantities 
with this mass. Thus we define T — T/m,^, rrir — mr/mi, and we set the renor- 
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malization scale /i equal to nii, (which, of course, does not affect any physics). We 
then consider the dimensionless effective potential 

v(6,6) = m;^(^tot(6,6)-v;ot(6 = ^,C3 = o)) 

= A{f)M^ + ^^{MY (2.52) 

OTT 

where A(f) = m;^ + fy6, B(f)^n- logf\ 

We can now examine the phase structure of V as a function of temperature. 
The first fact we can prove is that if ml > 0, then the ground state is always given 
by Al^ = 0. This is easiest to see by completing the square in V to write it as: 

V(a, 6) = -27r'^^ + ^{M' + 47r'^f (2.53) 

If To is the zero of B{T) (Tq = exp(n/2) ~ 0.34), then for T > Tq, the quantity 
g(T) = -A'K'^A{f)/B{f) is strictly greater than 2, while for f < % g{f) is 
negative. Using these facts, as well as the restriction that M.^ lies between and 
1 in plotting the parabola described by V as a function of M."^ proves our claim. 
Thus domain walls appear when the ^j's roll to their minima in the > case 
(if initial condidtions allow for this; see below). 

The < case is more interesting. Set a — —Gm^/TQ. In order that the high 
temperature approximation remain valid during our analysis, we restrict ourselves 
to vahics of T satisfying T >> 1. 

First consider the a < 1 case. Then for the values of T under consideration, 
A{f) > 0, B{f) < 0. Furthermore, the location of the peak of V is at Ml,^ ~ 
0.77(x — Q;)/logx, where x = (T/Tq)^. Now, T » 1 implies that we must only 
consider values of x larger than Tq"^ ~ 8.5. In this regime, TW^ax always larger 
than one, so that again, a plot of the parabola described by V as a function of Al^ 
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leads to the conclusion that the ground state is at = 0, just as in the > 
case. 

What happens in the a > 1 case? For a « Tq"^, we are essentially back 
to the situation described in the previous paragraph. However, suppose that a 
is comparable to or larger than Tq^. Then it is simple to see (using once more 
the parabola described by V as a function of Ai^, see fig. (2)) that the ground 
state switches from A^^ = 0to7W^ = lata critical temperature such that 
—Stt"^ A(Tc) / B (Tc) = 1. This is a second order phase transition which takes the 
domain walls present when = was the ground state and converts them into 
cosmic strings which are allowed when = 1 is the vacuum state. This provides 
a natural solution to the problem of having a long-lived network of domain walls 
in theories of late time phase transitions. We display some snapshots of the phase 
transition in fig. (3). 

We see then that, as stated in the introduction, the vacuum structure of this 
class of models is surprisingly rich in terms of the evolution of topological defects. 
We can also find large regions of parameter space for which the critical temperature 
of both the phase transition that generates the domain walls as well as of the one 
that converts them into strings is O{10m,j). Thus these models are well suited to 
being used for late time phase transition purposes. 

Let us now return to considering the observational data that provides the 
strongest motivation for a continued investigation of LTPT's. It has recently been 
pointed out by M. Schmidt et al. that comoving quasar space densities exhibit 
a strong peak at redshifts of 2 to 3||l|. In their chapter they plot quasar space 
density as a function of redshift (z) and also separately as a function of cosmic 
time. There is an unmistakable peak in the quasar space density in the region of 
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redshift 2 to 3. The authors further point out that the observed dechne in space 
density for space density for z > 3 is not a result of instrumentational difficulties in 
detecting distant quasars. On the contrary, the decline in observed quasar density 
is a real decline in the number density for 2; > 30. 

Thus, Schmidt et al. point out that quasar density peaks sharply at an epoch 
of about 2. 3 billion years after the Big Bang. The full width at half maximum of 
the peak is around 1. 4 billion years. Their discussion makes the final point that 
one needs to understand these time scales in terms of the formation and evolution 
of quasars. 

What we would like to argue here is that such a peaked distribution of objects 
may result from the formation of topological defects as the universe goes through 
a phase transition. Such topological defects could form the seeds around which 
quasars light up. 

The central power supply of quasars is believed to be gravitational in origin[§]. 
It is suggested in this chapter that this central power supply may be formed as 
a result of topological defect formation in a phase transition linked with massive 
neutrinos. Such a phase transition would happen at the right epoch if one believes 
the neutrino masses implied by the MSW solution to the solar neutrino problem. 
The production of black holes as a result of cosmological phase transitions has 
been discussed by various authors]^, ||, 0]. 

The idea of topological defects as seeds for structure formation has been around 
for some time. There are many excellent and recent reviews on the subject [^, 
T0| . The idea that topological defects formed after the decoupling of matter and 



radiation may play an important role in structure formation has also been discussed 
before by Hill, Schramm and Fry[|n[|(HSF). 
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In the Standard Cosmology the redshift range of 2 to 3 corresponds to a few 



meV in mass scalesp6|. Thus the critical temperature of the phase transition 



required to produce quasars is fixed at a few meV. 

As we have already pointed out in the context of our particle physics model for 
LTPT's we expect a phase transition with a critical temperature of a few times 
the relevant light neutrino mass. Let us therefore now turn to a discussion of 
what evidence we have for non-zero masses for neutrinos and what masses for light 
neutrinos are implied by these observations. 

At neutrino detectors around the world, fewer electron neutrinos are received 
from the sun than predicted by the Standard Solar Model. An explanation of the 
deficiency is offered by the MSW mechanism]^ which allows the v^. produced in 
solar nuclear reactions to change into z/^. This phenomenon of neutrino mixing 
requires massive neutrinos with the masses for the different generations different 
from each other Qj. 

The model we considered earlier was an extremely simple one. Although it 
had 2 families of light neutrinos, there was only one single light neutrino mass. 
As such this model was not compatible with the MSW effect. However it is fairly 
straightforward to modify our original model to make it compatible with the MSW 
effect as is shown in what follows. 

To ensure that it is not possible to choose the weak interaction eigenstates to 
coincide with the mass eigenstates we must require the 2 neutrino mass scales to 
be different. We can ensure neutrino mixing in our model by demanding that niai 
be such that m\\ ^ and mi2 = = m2i. In this case, the effective potential 
Kot(^i, ^3) = l/2{y{M\) + V{Ml)) with Vi^Mf) having the same functional form 
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a.sV{A4^){i = 1, 2) and A4j given by the following expression: 



M^ = ^(cos^2||e||+8sin2 2||e||) (2.54) 



Further, if mu « 777,22 then Vi;ot('Ci) ^3) = Vi-^Df^^ which is exactly half the 
finite temperature effective potential we discussed earlier except the neutrino mass 
scale is the heavier neutrino mass scale. Hence, the discussion on phase transitions 
and formation of topological defects we carried out earlier goes through exactly 
except that the critical temperature is determined by the mass scale of the heavier 
of the 2 neutrinos. 

In the complete picture of neutrino masses^], the neutrinos might have a mass 
hierarchy analogous to those of other fermions. Further, we expect that the mixing 
between the first and third generation might be particularly small . Thus we will 
only consider z/g-z/^ mixing. The data seems to imply a central value for the mass 



of the muon neutrino to be a few meV 27 



At this point, let us re-examine the constraint on the domain wall tension, a , 
placed by the measurements of the CMBR on small angular scales. Recall that the 
constraint is cr < 0.5MeV^ . An estimate of a in terms of the quantities rrii, and 
/ introduced in our model can be obtained0. (To make contact with the work 
of L. Widrow cited above please note that his A777^ = m^(z/^) and ?n = / in our 
notation. ) Thus, the constraint on a then implies that / < lO^^GeV. Our model 
is clearly an effective theory with / being some higher symmetry breaking scale on 
which it is tough to get an experimental handle. However, the constraint derived 
above is in fact natural in the context of the see-saw model of neutrino masses 
embedded in Grand Unified Theories as discussed by Mohapatra and Parida| 
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and also by Deshpande, Keith and Palp9 



There are a few points on which we'd hke to add some further comments. The 
first has to do with black hole formation as a result of LTPT's. The second with 
the possible role of a LTPT linked with the mass of v-r in the formation of LSS (on 
scales ~ lOOMpc). Finally, a brief comment on the post COBE status of LTPTs 
is made. 

A number of different groups have studied different mechanisms for black hole 
formation as a result of cosmological phase transitions. Black hole formation as a 
result of bubble wall collisions has been suggested by Hawking et al[0. Collapse of 
trapped false vacuum domains to produce black holes has been studied by Kodama 
et al.[0 and Hsu[^]. The collapse of closed domain walls to form black holes has 
been studied by Ipser and Sikivie[^] and also by Widrowf^. In fact, Widrow has 
suggested that the collapse of closed domain walls is likely to produce black holes 
as a result of LTPTs. An estimate of the mass of black holes produced in LTPTs 
is Mbh ~ 10^ Mgofar ( iqs A/pc )^ ( w^GeV ) ' "^^erc R is the radius of the closed domain 
wall and we have taken m(i/^) to be 3meV to obtain our estimate. 

There are of course two definite uncertainities in getting a number out of the 
above expression. First, there will obviously be a range of possible sizes of closed 
domain wall bubbles that will be formed in the transient period of the phase 
transition. A natural upper limit to the size of the bubble will be the horizon 
size at the epoch of the phase transitionThus this implies that R < lO^Mpc. The 
horizon size at the onset of the phase transition is smaller. Further, it should be 
noted that sub-horizon sized bubbles which will be formed in the transient period of 
the phase transition may also play an important role. Secondly, as already pointed 
out our model is clearly an effective theory with / being some higher symmetry 
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breaking scale on which it is tough to get an experimental handle. 

A more detailed analysis of the masses of BH produced in LTPTs as well as the 
number densities of these BH will be the subject of a later work. For now, the fact 
that it is possible to produce BH which may act as the central engine of quasars is 
pointed out. One can check that the mechanism of collapse of domain wall bubbles 
can in fact give roughly the needed number of black holes to power quasars. The 
Hubble radius at the onset of the phase transition at z ~ 5 , is a few hundred Mpc. 
One expects a distribution of bubble sizes with the horizon as an upper limit to 
be formed during a cosmological phase transition. In fact, Turner, Weinberg and 
Widrow[^ have carried out a detailed study of the distribution of bubble sizes 
resulting from cosmological phase transitions. They report that typical bubble sizes 
in a successful phase transition range from 0.01 to 1 times the Hubble radius at the 
epoch of the phase transition and depends only very weakly on the energy scale of 
the phase transition. Thus, the most massive black holes formed from the phase 
transition being studied here will have an abundance ~ 10~^ to lQ~^Mpc~^ with a 
higher abundance of less massive black holes. These numbers, in fact match very 
well with the observational number densities of quasars as discussed by Warren and 
Hewett[p4|, and by Boyle et al. and Irwin et al|35]. The point is that the epoch 



of the phase transition linked with massive neutrinos is about right to explain the 
observed peak in the quasar distribution. Further, order of magnitude estimates 
of the masses of black holes produced and their abundances are consistent with 
observations within the uncertainities discussed above. A detailed analysis of the 
efficiency of black hole formation and therefore an accurate number for the space 
density of black holes produced in this model will be carried out in a later work. 
The other point that needs to be commented on is the role of a massive Vr in 
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the scheme of things presented here. The mass of Vr is less well determined but 
one would expect an earlier phase transition with a Tc of a few times m{vr). In 
fact, such a phase transition may well be responsible for the LSS seen on the scale 
of ~ lOOMpc. There already exists a great deal of discussion on this subject. 

However, black hole formation in this earlier phase transition would be more 
difficult to observe. In fact, it is the black holes that are formed in the most recent 
phase transition that will have the greatest observable consequences. 

The anistropy of the CMBR on large angular scales is more closely linked to 
the formation of LSS (on scales ~ 100 Mpc). One can relate the power spectrum 
of density fluctuations to the gravitational potential power spectrum responsible 
for distortions of the CMBR. This link is thoroughly discussed by Jaffe, Stebbins 
and Frieman (JSF)|3^. Since considerable processing of the power spectrum must 
take place in the process of black hole formation a relationship between quasar 
distribution and the distortions of the CMBR on large angular scales is consider- 
ably more difficult to establish. The viability of LTPTs in the post-COBE era has 
recently been discussed by Schramm and LuofS^. They point out that LTPTs are 
still a viable model for the formation of LSS. 

To place things in perspective one should keep in mind that different mecha- 
nisms may play important roles in structure formation at different length scales. 
This point has been emphasised by Carr[^who has given a comprehensive discuss- 
sion on the origin of cosmo logical density fluctuations. In fact, JSF also conclude 
their discussion by pointing out that the final power spectrum is most likely due 
to a combination of primordial and late time effects. 
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2.4 Conclusion 



This chapter provides a detailed apphcation of the imaginary time formahsm of 
finite temperature field theory. In addition to providing an illustration of many 
useful techniques this is also of current cosmological interest because of a number 
of recent observations and measurements that it has the promise of tying up. Thus, 
the MSW solution to the solar neutrino problem seems to imply a muon neutrino 
mass of a few meV. This in turn would lead to a phase transition in the PNGB 
fields associated with massive neutrinos with a critical temperature of several meV. 
This phase transition happens at the correct epoch in the evolution of the universe 
to provide a possible explanation of the peak in quasar space density at redshifts 
of 2 to 3. This work is clearly only a first step in bringing these ideas together. 
There are clearly some issues which need to be addressed in greater depth and 
explored in fuller detail. Thus, a more detailed and accurate analysis of the black 
hole formation efficiency and the masses and number densities of BH's needs to be 
carried out. 

To obtain more precise estimates of the black hole number densities and masses 
one needs to make a more detailed analysis of domain formation and growth in 
FRW cosmologies. Domain formation and growth are among the most impor- 
tant non-equilibrium phenomena associated with phase transitions. The real time 
formulation of finite-temperature field theory is the most logical and well-suited 
formalism to study such phenomena. 

Boyanovsky, Lee and Singh have studied domain formation and growth and 
discussed it at length using the real time time formalism[0. This study was 
carried out in Minkowski space. We need to extend this work to FRW cosmologies. 

60 



In fact, non-equilibrium phenomena such as particle creation, entropy growth and 
dissipation have already been studied for FRW cosmologies by Boyanovsky, de 
Vega and Holman|[40||. We now need to examine domain formation and growth in 
this formalism. 

This number density and mass distribution of BH's then needs to be compared 
to the observational data on quasars. Further, a more detailed analysis of the 
angular dependence of the CMBR anisotropy produced needs to be carried out 
and compared to observations. 

Another interesting cosmological implication of the calculation outlined above 
is the potential it has of resolving the discrepancy between the recently measured 
values of the Hubble constant and the age of the universe. Increasing improvements 
in the independent determinations of the Hubble constant and the age of the 
universe now seem to indicate that we need a small non-vanishing cosmological 
constant to make the two independent observations consistent with each other. 
The cosmological constant can be physically interpreted as due to the vacuum 
energy of quantized fields. To make the cosmological observations consistent with 
each other we would need a vacuum energy density, p„ ~ (lO^^el^)'* today ( in the 
cosmological units h = c = k = l). In a recent paper we have argued that such 
a vacuum energy density is natural in the context of phase transitions linked to 
massive neutrinos. In fact, the neutrino masses required to solve the cosmological 
constant problem are consistent with those required to solve the solar neutrino 
problem by the MSW mechanism. We will return to a discussion of this subject in 
the last chapter. My hope is that this work will stimulate further thoughts about 
these issues. 
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Figures for Chapter 2: 

Figure 1: V as a function for Ad"^ for temperatures just above, at and just below 
the phase transition temperature. 

Figure 2: V(^i,^3) for the temperatures used in fig. (2). 



67 



Chapter 3 



REAL TIME APPLICATION: 
DISSIPATION VIA PARTICLE 
PRODUCTION IN SCALAR 
FIELD THEORIES 

3.1 Introduction 

In this chapter we will study the dynamics of dissipation in scalar field theories 
using the real time formalism. Dissipation and reheating are among the most 
important non-equilibrium phenomena that accompany any phase transition. In 
cosmology these process are particularly important in the context of infiationJI], 0. 
Recall that at the end of new or chaotic infiation^, defined by when the slow-roll 
conditions for the so-called infiaton field0, ^, (f){t), fail to obtain, the infiaton 
begins to oscillate about its true ground state. 

The dissipation of the coherent oscillations of the field due to particle produc- 
tion transfers the enormous vacuum energy of the fields prior to the phase transition 
into the energy of the particles. The subsequent thermalisation of these particles 
allows the infiationary cosmology to smoothly merge into the standard hot big 
bang cosmology. Futher, the enormous creation of entropy which accompanies the 
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particle production helps inflation solve the problems of standard cosmology. 

In all realistic models of inflation, the inflaton is coupled to both fermionic 
and bosonic lighter flelds. As the inflaton oscillates about the minimum of its 
potential, it decays to the lighter flelds it is coupled to. This process is usually 
modeled by taking the evolution equation of the inflaton cf) and inserting by hand 
a term which looks like r0|^, |^. It is important to state that the Tcj) term is not 
derivable starting from the fundamental Lagrangian or Hamiltonian of the system 
but has to be put in by hand just to model the dissipative dynamics. A more 
careful treatment is both desirable and required and this has been pointed also by 
Linde et al and Brandenberger et al]^ 

This chapter will lay out the formalism that allows us to study the detailed 
dynamics of dissipation via particle production starting from the exact microscopic 
Hamiltonian of the system. 

Though we have studied the time evolution of a scalar fleld coupled to fermions 
as well as another scalar fleld, in this chapter we will study the dynamics of a self- 
interacting scalar fleld. This is a natural flrst step. It will allow us to lay out the 
formalism in the simplest possible setting and allow us to gain insight into some 
of the essential physics in this simple setting. 

One of the issues which frequently arises in trying to understand dissipation 
from a microscopic perspective is that dissipation appears to be an irreveresible 
process whereas the underlying microscopic evolution equation are time reversible. 
In this context, it is important to realize the importance of two important factors 
that alow us to obtain dissipative dynamics from time reversal invariant evolution 
equations. These two factors are non-equilibrium initial conditions and coarse 
graining. Thus, even though the evolution equations are time reversal invariant, the 
non-equilibrium initial conditions make the behaviour of the system irreversible. 
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Further, coarse graining implies that even though total energy is conserved there 
is a transfer of energy from the modes of interest into the remaining modes which 
now play the role of an environment to the field modes or interest are coupled to. 

Here is the plan for this chapter. For a self-interacting scalar field theory we 
study the time evolution of system by having the expectation value of the field away 
from the minimum of its potential. This provides the necessary non-equilibrium 
initial condition. The expectation value of the field starts oscillating about the 
minimum. These oscillations of the expectation value start damping out as the 
fluctuations start building up. We also study particle production taking place in 
the system. Thus, we arrive at the conclusion that there is a precise correlation 
between the damping out of the oscillations of the zero mode with the building up 
of the fluctuations and particle production. 

We study two different self-interacting scalar field theories. One possessing a 
discrete symmetry and another possessing a continuous symmetry. For both these 
systems we investigate what happens when the symmetry is unbroken and then 
when the symmetry is broken. We will notice that the most dramatic damping 
takes place in the broken continuous symmetry case because of the availability of 
the Goldstone channel for dissipation. 

We start off by setting up the formalism to study the time evolution. The 
time evolution of the system is determined by the quantum Liouville equation. 
We are interested in the time evolution of the expectation value of the field. (This 
time evolution equation is in fact a functional form of the Ehrenfest Theorem.) 
The field can be decomposed into the expectation value and fluctuations about 
the expectation value. As a flrst step we can try to do a loop expansion. How- 
ever, we will notice that very soon the one-loop contribution becomes of the same 
order as the zeroth order contribution and hence one cannot trust the perturba- 
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tive expansion beyond this time. Thus, one is driven to consider ways to extract 
non-pcrturbativc information. We next use the self-consistent Hartree Approxima- 
tion to extract information of non-perturbation dynamics. We will then introduce 
the formalism to study particle production. The evolution equation we arrive at 
are coupled integro-differential equation which we solve numerically. Plots of the 
time evolution of the expectation value of the field, the growth of fiuctuations and 
particle production will allow us to gain insight into the physics of this problems. 
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3.2 Non-equilibrium field theory and equations 
of motion 



The generalization of statistical mechanics techniques to the description of non- 



equilibrium processes in quantum field theory has been available for a long time[|ri 
[T2I , [13| , [T^ , |T5[] but somehow has not yet been accepted as an integral part of the 
available tools to study field theory in extreme environments. We thus begin 
by presenting a somewhat pedagogical introduction to the subject for the non- 
practitioner. The discussion that follows is based on understanding gained through 



collaborations with Boyanovosky, de Vega, Holman and Lee|T^, |T^, |T8|, |T9|, pO 

The non-equilibrium description of a system is determined by the time evolution 
of the density matrix that describes it. This time evolution (in the Schroedinger 
picture) is determined by the quantum Liouville equation. A non-equilibrium 
situation arises whenever the Hamiltonian does not commute with the density 
matrix. Here we allow for an explicitly time dependent Hamiltonian, which might 
be the case if the system is in an external time dependent background, for example. 
Our starting point is the Liouville equation for the density matrix of the system: 

,n^ = [H,^,,,pit)], (3.1) 

where -ffevoi is the Hamiltonian of the system that drives the non- equilibrium 
evolution of the system. 

Next we use the density matrix to extract the order parameter from the field 
in the Schroedinger picture: 

0W = ^/ d'xTTipmm, (3.2) 

with Q the spatial volume (later taken to infinity), and $(x) the field in the 
Schroedinger picture. Using the Liouville equation together with the Hamiltonian: 
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H = /rf^x|^^ + ^(V$)^ + y($)| (3.3) 

and the standard equal time commutation relations between a field and its canon- 
ically conjugate momentum, we find the equations: 



^ = i/rfMn(x,t)) = i/rf3xTr[p(t)n(f)]=7r(t) (3.4) 
From these equations we can find the equation of motion for the order parameter 

m 

We expand the field operator as $(x) = (f){t) + ilj{x, t), insert this expansion into 
equation ( |3.6| ) and keep only the quadratic terms in the fluctuation field %l){x,t). 
Doing this yields the equation: 



''^^l*) + V'm) + / + . . . = (3.7) 



Here the primes stand for derivatives with respect to cj). 

At this point we need to compute and keep track of the fluctuations in the field. 
The fluctuations {ilj^{x,t)) is the two-point correlation function. In the language 
of field theory this is the two-point equal time Green's function. This Green's 
function can be evaluated in terms of the homogeneous solutions of the quadratic 
form for the fluctuations. For concreteness let us consider a self- interacting scalar 
field with the potential, 

F($) = lm2$2 + A$4 ^3_g) 
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In this case the homogeneous solutions of the quadratic form for the fluctuations 
are detemined by 







with 



(3.9) 
(3.10) 

(3.11) 



The boundary conditions (|3.1CI|) correspond to positive and negative U~ fre- 
quency modes for t < (the Wronskian of these solutions is 2iujl). Notice that 
U^{t) = [U;l(t)]*. In terms of these mode functions we obtain the two-point equal 
time Green's function to be 



G{t,t) 



(3.12) 



Thus to order h we find the following equations 
0(t) + + -03(,) + _^(,) J 



^ + P + m^ + ^0^(t) 



(3.13) 

(3.14) 
(3.15) 



where we have restored the Ti to make the quantum corrections explicit. This set of 
equations clearly shows how the expectation value (coarse grained variable) "trans- 
fers energy" to the mode functions via a time dependent frequency, which then in 
turn modify the equations of motion for the expectation value. The equation for 



the mode functions, (|3.14|) may be solved in a perturbative expansion in terms of 
A0^(t) involving the retarded Green's function. 

Before attempting a numerical solution of the above equations it is important 
to understand the renormalization aspects. For this we need the large k behaviour 
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of the mode functions which is obtained via a WKB expansion as in references 



2^ , |T9[ and to which the reader is referred to for details. We obtain 

d'k |f/+(t)p 



(27r 



)3 2a;0 



A2 
Stt 



in 



m2 + ^02(t) 



In 


A- 






+ 









+ finite 



(3.16) 



where A is an ultraviolet cutoff and k an arbitrary renormalization scale. From the 
above expression it is clear how the mass and coupling constant are renormalized. 
It proves more convenient to subtract the one-loop contribution at t = and absorb 
a finite renormalization in the mass, finally obtaining the renormalized equation 
of motion 

)u^{t)\'-i\ 



OTT^ Jo 



k'dk 



00 







+ ^0W(0'W-0'(O))ln[A/K] 



(3.17) 



In the equations for the mode functions the mass and coupling may be replaced 
by the renormalized quantities to this order. One would be tempted to pursue 
a numerical solution of these coupled equations. However doing so would not 
be consistent, since these equations were obtained only to order h and a naive 
numerical solution will produce higher powers of Ti that are not justified. 

Within the spirit of the loop expansion we must be consistent and only keep 
terms of order h. First we introduce dimensionless variables 



X 



6m^ 



k 

(r) ; T = niRt ; q = 

rriR 



9 



87r2 



(3.18) 



and expand the field in terms of g as 

Vi^) = Vci{r) + gvi{r) + ■ ■ ■ 
Now the equations of motion consistent up to 0{h) become 
Vci{r) + r]4T) + r]l{T) = 



(3.19) 



(3.20) 
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Jo 



m{r)\^-l 



[g2 + l + 3r/2(0)] 



+ 



Ivciir) {vUr) - vliO)) H^/mn] = (3.21) 
The solution to equation ( |3.20|) is an elliptic function. The equations for the 



mode functions become 



+ + I + ^Vliir) 



Kir) = (3.22) 



with the boundary conditions as in eq. ( |3.15|) in terms of the dimensionless frequen- 



cies and where for simplicity, we have chosen the renormalization scale k = m^. 
The chosen boundary conditions //(O) = r/o ; ^(0) = can be implemented as 

r/e/(0) = r/o ; 7)^(0) = ; r/i(0) = ; 7)1(0) = (3.23) 

In fig. (2) we show rjiij) with the above boundary conditions with 770 = 1 and 
g = 0.1. A cutoff A/mji = 100 was chosen but no cut-off sensitivity was detected 
by varying the cutoff by a factor 3. Notice that the amplitude grows as a function 
of time. This is the behavior shown in fig. (2). 

Failure of perturbation theory to describe dissipation 

This section has been devoted to a perturbative analysis of the "dissipative as- 
pects" of the equation of motion for the scalar field. Perturbation theory has been 
carried out as an expansion up to 0{h), both in the broken and the unbroken sym- 
metry case. The potential for the unbroken and broken symmetry cases is given 
in figures 3 and 4 respectively. In both cases we found that the amplitude of the 
quantum corrections grow as a function of time and that the long time behavior 
cannot be captured in perturbation theory. This failure of perturbation theory to 
describe dissipation is clearly understood from a very elementary but yet illumi- 
nating example: the damped harmonic oscillator. Consider a damped harmonic 
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oscillator 

q + rq + q = (3.24) 

with r ~ 0{X) where A is a small perturbative coupling. The above equation can 
be solved exactly to yield the solution, 

r r 

q{t) = e-2* cos[cu(r2)t] ^ cos(t) - -t cos(t) + ^(r^). (3.25) 

We see that there does exist a perturbative expansion in F, but in order to find 
appreciable damping, we must wait a time ~ 0{l/r) at which perturbation theory 
becomes unreliable. 

In order to properly describe dissipation and damping one must resum the 
perturbative expansion. Another hint that points to a resummation of the per- 
turbative series is provided by the set of equations eqs. (|3.2CI| - |3.22|) . In eq.( p.20|) . 



the classical solution is a periodic function of time of constant amplitude, since 
the classical equation has a conserved energy. As a consequence, the "potential" 
in the equation for the modes (eq.( |3.22]) ) is a periodic function of time with con- 
stant amplitude. Thus although the fiuctuations react back on the coarse grained 
field, only the classically conserved part of the motion of the coarse grained field 
enters in the evolution equations of the mode functions. This is a result of being 
consistent with the loop expansion, but clearly this approximation is not energy 
conserving. 

As we will point out in the next section, in an energy conserving scheme the 
fiuctuations and amplitudes will grow up to a maximum value and then will always 
remain bounded at all times. 

Thus in summary for this section, we draw the conclusion that perturbation 
theory is not sufficient (without major ad-hoc assumptions) to capture the physics 
of dissipation and damping in real time. A resummation scheme is needed that 
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effectively sums up the whole (or partial) perturbative series in a consistent and/or 
controlled manner, and which provides a reliable estimate for the long-time behav- 
ior. The next section is devoted to the analytical and numerical study of some of 
these schemes. 

3.3 Non-perturbative schemes I: Hartree approx- 
imation 

Motivated by the failure of the loop expansion, we now proceed to consider the 
equations of motion in some non-perturbative schemes. First we study a single 
scalar model in the time dependent Hartree approximation. After this, we study 
an 0(N) scalar theory in the large N limit. This last case allows us to study the 
effect of Goldstone bosons on the time evolution of the order parameter. 

In a single scalar model with the interactions described by the potential of 
eq. (|3.8|) , the Hartree approximation is implemented as follows. We again decom- 
pose the fields as $ = (p + ip. The Hartree approximation is implemented by 
factorising all powers of the field higher than two and expressing them in terms of 
products of two point correlation functions and upto quadratic terms in the field 
and demanding that a self-consistency condition be satisfied. Thus, the Hartree 
approximation is obtained by assuming the factorization 

^ 3(V''(x,t))^(f,t) 
^^{x,t) 6{ip^{x,t))ip^{x,t) -3{'4j^{x,t))^ (3.26) 

Translational invariance shows that {tp'^{x,t)) can only be a function of time. The 
expectation value will be determined within a self-consistent approximation. 

In this approximation the resulting Hamiltonian is quadratic, with a linear term 
in ip: 
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HMt)^J^^.{^.^^^mt).^^j. (3.27) 

Here 11 is the canonical momenta conjugate to iIj{x) and is given by, 



v«(t) = n0) + ^A0(v^2) 



(3.28) 



It is recognized as the derivative of the Hartree "effective potential" ||2^, ^ 
with respect to (it is the derivative of the non-gradient terms of the effective 



action ll, |T8|, |25[). 
Also, 



(3.29) 



With this Hartree factorization we can use the non-equilibrium formalism outlined 
in the previous section to determine the time evolution. The resulting equations 
are obtained to be: 



2 



mt)) 



-tGk{t,t)] 



d?k 




(27r)3 


2cu,(0) 



(27r)3 

t/+(t)=0 ; ul{t) = k' + M\t) 



(3.30) 

(3.31) 
(3.32) 



The initial conditions for the mode functions are 

[/+(0) = 1 ; f7+(0) = -z^fc(O) 



(3.33) 



It is clear that the Hartree approximation makes the Lagrangian density quadratic 
at the expense of a self-consistent condition. In the time independent case, this 
approximation sums up all the "daisy" (or "cactus") diagrams and leads to a 
self-consistent gap equation. 
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At this stage, we must point out that the Hartree approximation is uncontrolled 
in this single scalar theory. This approximation does, however, become exact in 
the N —>■ oo limit of the 0{N) model which we will discuss in the next section. 

3.3.1 Renormalization 

We now analyze the renormalization aspects within the Hartree approximation. To 
study the renormalization we need to understand the divergences in this integral 







(27r)3 





(3.34) 



The divergences will be determined from the large-k behavior of the mode func- 
tions which obey the differential equations obtained from ( p.32| ) with the initial 
conditions (|3.33|) . By a WKB-type analysis (see|2^, |l9l for a detailed description), 
in the k ^ oo limit, we find 



2cufc(0) 



1 1 



(3.35) 



Inserting these results in the integral, it is straightforward to find the divergent 
terms and we find 



d?k 


utit) " 


(27r)3 


2^fc(0) 



57r^ V K 



+ finite 
(3.36) 

where A is an upper momentum cutoff and k a renormalization scale. 

Now, we are in position to specify the renormalization prescription within the 
Hartree approximation. In this approximation, there are no interactions, since 
the Lagrangian density is quadratic. The nonlinearities are encoded in the self- 
consistency conditions. Because of this, there are no counterterms with which to 
cancel the divergence and the differential equation for the mode functions must be 
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finite. Therefore, it leads to tlie following renormalization prescription: 



ml + + =mi+ + ^mt))R (3.37) 

where the subscripts B, R refer to the bare and renormalized quantities respectively 
and {i''^{t))B is read from eg. ( p.36| ): 



-^A ^ In 



X 



X 



^ /„/.2/ 



+ finite (3.38) 



Using this renormalization prescription eq.( 3.37|) , we obtain 



2 . -^B .2 2 



167r2 



A 



B 



(3.39) 
(3.40) 



and 



X 



/ 







(27r)3 


2c^fc(0) 



:A^ 



57r 



It is clear that there is no wavefunction renormalization. This is a consequence of 
the approximation invoked. There is, in fact, no wavefunction renormalization in 
either one-loop or Hartree approximation for a scalar field theory in three spatial 
dimensions. 

With an eye towards the numerical analysis, it is more convenient to write 



R 



(3.41) 



and perform a subtraction at time t = absorbing {iIj^{0))b into a further finite 
renormalization of the mass term (m| + (^^(0))/?, = M|). 
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The renormalized equations that we will solve finally become 



1 



(3.42) 

U^(t) = (3.43) 



and 



1 - ^ In ( A 



X 



r d"! 
J (2^ 



(27r)3 2a;fe(0) IGtt^ 
with the initial conditons for U^{t): 



in(^) [m-m) 



t/+(0) = 1 ; f7+(0) = -z^fc(O) ; cu,(0) = s^k^ + Ml+ ^4>\Q) (3.44) 

It is worth noticing that there is a weak cutoff dependence on the renormalized 
equations of motion of the order parameter and the mode functions. This is a con- 
sequence of the well known "triviality" problem of the scalar quartic interaction in 
four space-time dimensions. This has the consequence that for a fixed renormal- 
ized coupling the cutoff must be kept fixed and finite. The presence of the Landau 
pole prevents taking the limit of the ultraviolet cutoff to infinity while keeping the 
renormalized coupling fixed. 

This theory is sensible only as a low-energy cutoff effective theory. We then 
must be careful that for a fixed value of Ar, the cutoff must be such that the theory 
never crosses the Landau pole. Thus from a numerical perspective there will always 
be a cutoff sensitivity in the theory. However, for small coupling we expect the 
cutoff dependence to be rather weak (this will be confirmed numerically) provided 
the cutoff is far away from the Landau pole. 
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3.3.2 Particle Production 



Before we engage ourselves in a numerical integration of the above equations of 
motion we want to address the issue of particle production since it is of great 
importance for the understanding of dissipative processes. 

In what follows, we consider particle production due to the time varying effec- 
tive mass Ai'^it) in eq.( p.29| ) of the quantum field ip for the single scalar model. 

In a time dependent background the concept of particle is ambiguous and it 
must be defined with respect to some particular state. Let us consider the Heisen- 
berg fields at t = written as 



^(x, 0) = ^ E / V , («^-(0) + alM 



n^(f,0) = («^(0) - alM) e^-' (3.45) 



with co'fc(O) as in eq.( p.32| ) and the spatial volume. The Hamiltonian is diagonal- 
ized at t = by these creation and destruction operators: 



i7(0) = E^/.(0) 

k 



4(0)afc(0) + i 



(3.46) 



Thus we define the Hartree-Fock states at t = as the vacuum annihilated by 
afc(O) together with the tower of excitations obtained by applying polynomials in 
aj,(0) to this vacuum state. The Hartree-Fock vacuum state at t = is chosen as 
the reference state. As time passes, particles (as defined with respect to this state) 



will be produced as a result of parametric amplification p6|, p7| . We should mention 



that our definition differs from that of other authors pT], p7[] in that we chose the 
state at time t = rather than using the adiabatic modes (that diagonalize the 
instantaneous Hamiltonian). 
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We define the number density of particles as a function of time as 



d^k Tr [4(0)afc(0)p(t) 



(27r)3 
where by definition 



Trp(O) 



d^k Tr 4(t)afc(t)p(0) 



(2vr) 



Trp(O) 



(3.47) 



4W = U-\t,0)al{0)U{t,0) ; ak{t) = U-\t,0)ak{0)U{t,0) 



(3.48) 



are the time-evolved operators in the Heisenberg picture. In terms of these time 
evolved operators, we may write: 



ip{x, t) 

n^(f, t) 



—I 



(3.49) 



On the other hand, we now expand the Heisenberg fields at time t in the 
following orthonormal basis 



ip{x, t) 

n^(f,t) 



E 

k 

E 



1 



2cufc(0) 



~akUt{t) + ~alj,U^{t))e 



Ak-x 



~akU^{t) + dl,U,{t))e 



Ak-x 



(3.50) 



k 

where the mode functions U'^it) ; Uj^(t) = (f/^(t))* are the Hartree-Fock mode 
functions obeying eqs.( p.43| - p.44| ) together with the self consistency condition. 

Thus Sfc, dl are the annihilation and creation operators of Hartree-Fock states, 
and the Heisenberg field ipi^^ t) is a solution of the Heisenberg equations of motion 
in the Hartree approximation. Therefore the dl and dk do not depend on time and 
are identical to a|.(0) and 0^(0) respectively ( we can check this by evaluating the 
expansion in eq. (|3.50|) at t = together with the initial conditions on U'^{t) in 
eq. (|3.33| )). Using the Wronskian properties of the function Uj^it), we see that the 
dl and d^ satisfy the usual canonical commutation relations. The reason for the 
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choice of the vacuum state at t = now becomes clear; this is the initial time at 
which the boundary conditions on the modes are determined. The mode functions 
f/^(t) ; U^{t) are then identified with positive and negative frequency modes at 
the initial time. 

By comparing the expansion in eg. ( p.45| ) evaluated at time t with that in 
eq.( p.50| ), we find that the creation and annihilation operator at time t can be 
related to those at time t = via a Bogoliubov transformatiom: 



ak{t) = J^+,k{t)ak{0) + ^_,fe(t)alfc(0) 
The J-'±{t) can be read off in terms of the mode functions U^{t) 



(3.51) 



\ Kit) ' 



1 + 







2 







+ 



\T+,k{tt - l-^-,A.(t)|' = l 



(3.52) 



At any time t the expectation value of the number operator for the quanta of 
ip in each k-mode is given by 

'al{t)ak{t)p{0) 



After some algebra, we find 



Tr 



Trp{0) 



(3.53) 



A4(t) = (2 - 1) A4(0) + (|^+,fc(t)f 



(3.54) 



This result exhibits the contributions from "spontaneous" ( proportional to the 
initial particle occupations) and "induced" ( independent of it) particle production. 
Since we are analyzing the zero temperature case with A/fc(0) = only the induced 
contribution results. 
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3.3.3 Numerical Analysis 
Unbroken Symmetry Case 

In order to perforin a numerical analysis it is necessary to introduce dimensionless 
quantities and it becomes convenient to chose the renormalization point k = Mr- 
Thus we define 



r]{t) = m\ 



A 



R 



2M% 



; Q 



M, 



k , , Xr 



R 



47r2 
Ml 



'Wr - mO))n) 



and finally, the equations of motion become 



d 



r] + r) + 



1 - 



d 



U+{t) = 



T)^ + gr]T.{T) = 



[/+(0) = 1 ; ^KiO) = -t^q^ + 1 + r^2(0) 



(3.55) 



(3.56) 
(3.57) 



E(r) 



/ Q dq 
Jo 



X 



^ +lln(^]{rj\r)-rj\0))h.58) 



V + l + r/2(0) 2 \Mr 

In terms of the dimensionless quantities we obtain the number of particles 
within a correlation volume N{t) = J\f{t)/M^ 



1 /-A/Mh 



Kir) 



+ 



Kir) 



(3.59) 



^^2 + 1 + ^2(0) 

Figures (5.a,b,c) show /^(t), S(t) and A^(t) in the Hartree approximation, 
for g — 0.1, r}{0) — 1.0 and A/Mr — 100; we did not detect an appreciable 

cutoff dependence by varying the cutofi' between 50 and 200. Clearly there is 
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no appreciable damping in r){T). In fact it can be seen that the period of the 
oscillation is very close to 27r, which is the period of the classical solution of the 
linear theory. This is understood because the coefficient of the cubic term is very 
small and g^{T) is neghgible. Particle production is also neghgible. This situation 
should be contrasted with that shown in figures (6.a-c) and (7.a-c) in which there 
is dissipation and damping in the evolution of r}{T) for 7^(0) — 4, 5 respectively 
and the same values for g and the cutoff. There are several noteworthy features 
that can be deduced from these figures. First the fiuctuations become very large, 
such that ^^(r) becomes 0{1). Second, figures (6. a) and (7. a) clearly show that 
initially, channels are open and energy is transferred away from the q — mode 
of the field. Eventually however, these channels shut off, and the dynamics of 
the expectation value settles into an oscillatory motion. The time scale for the 
shutting off of the dissipative behavior decreases as 77(0) increases; it is about 25 
for r]{0) = 4, and about 18 for 77(0) = 5. This time scale is correlated with the 
time scale in which particles are produced by parametric amplification and the 
quantum fluctuations begin to plateau (flgures (6.b,c) and (7.b,c)). Clearly the 
dissipative mechanism which damps the motion of the expectation value is particle 
production. Furthermore, the long time dynamics for the expectation value does 
not correspond to exponential damping. In fact, we did not flnd any appreciable 
damping for r > 70 in these cases. It is illuminating to compare this situation 
with that of a smaller coupling depicted in figures (8.a-c) for 77(0) = 5 g = 0.05. 
Clearly the time scale for damping is much larger and there is still appreciable 
damping at r fa 100. Notice also that gT.{r) ^ 2 and that particles are being 
produced even at long times and this again correlates with the evidence that the 
field expectation value shows damped motion at long times, clearly showing that 
the numerical analysis has not reached the asymptotic regime. 
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The fundamental question to be raised at this point is: what is the origin of the 
damping in the evolution of the field expectation value? Clearly this is a collision- 
less process as collisions are not taken into account in the Hartree approximation 
(although the one-loop diagram that enters in the two particle collision amplitude 
with the two-particle cut is contained in the Hartree approximation and is respon- 
sible for thresholds to particle production). The physical mechanism is very similar 
to that of Landau damping in the coUisionless Vlasov equation for plasmas and 



also found in the study of strong electric fields in reference 1 10]. In the case un- 
der consideration, energy is transfered from the expectation value to the quantum 
fluctuations which back-react on the evolution of the field expectation value but 
out of phase. This phase difference between the oscillations of 77^ (t) and those of 
S(r) can be clearly seen to be vr in figures (6.a,b), (7.a,b) since the maxima of 
rj^ir) occur at the same times as the minima of S(r) and vice versa. 

This is an important point learned from our analysis and that is not a priori 
taken into account in the usual arguments for dissipation via collisions. The process 
of thermalization, however, will necessarily involve collisions and cannot be studied 
within the schemes addressed in this paper. 

Broken Symmetry 

The broken symmetry case is obtained by writing M|. = — < and using the 
scale fiji instead of M/j to define the dimensionless quantities as in eq.( |3.55D and 
the renormalization scale. The equations of motion in this case become 

-2\ 1 



1 



3/ 1-flnf^ 
2 Vmj? 



^ + q'-l + r]\T)+gJ:iT) 
d 

d^' 



T]' + gr]i:ir) = 



f/+(r) = (3.60) 



t/+(0) = 1 ; -[/+(0) = -t^q^ + 1 + 7^2(0) (3.61) 



with S(r) given in eq.( |3.58| ) but with Mr replaced by /ij?. The broken symmetry 
case is more subtle because of the possibility of unstable modes for initial conditions 
in which //(O) << 1. We have kept the boundary conditions of eq.( ^.61| ) the same 



as in eq.( |3.57| ). This corresponds to preparing an initial state as a Gaussian state 
centered at 77(0) with real and positive covariance (width) and letting it evolve for 



t > in the broken symmetry potential |2^, The number of particles produced 



within a correlation volume (now (ij^) is given by eq. (|3.59|) with Mr replaced by 

Figures (9.a-c) depict the dynamics for a broken symmetry case in which //(O) = 
10"^ i.e. very close to the top of the potential hill. Notice that as the field 
expectation value rolls down the hill the unstable modes make the fluctuation 
grow dramatically until about r ^ 50 at which point 5'S(r) ^ 1. At this point, the 
unstable growth of fluctuations shuts off and the field begins damped oscillatory 
motion around a mean value of about ~ 1.2. This point is a minimum of the 
effective action. The damping of these oscillations is very similar to the damping 
around the origin in the unbroken case. Most of the particle production and 
the largest quantum fluctuations occur when the field expectation value is rolling 
down the region for which there are unstable frequencies for the mode functions 
(see eq. (|3.60| )). This behavior is similar to that found previously by some of these 



authors [^ . 

3.4 Non-perturbative schemes II: Large N limit 
in the 0{N) Model 

Although the Hartree approximation offers a non-perturbative resummation of 
select terms, it is not a consistent approximation because there is no a priori small 
parameter that defines the approximation. Furthermore we want to study the 
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effects of dissipation by Goldstone bosons in a non perturbative but controlled 
expansion. 

In this section, we consider the 0{N) model in the large N limit. The large 
N limit has been used in studies of non-equilibrium dynamics PB| Q 
and it provides a very powerful tool for studying non-equilibrium dynamics non- 
perturbatively in a consistent manner. The Lagrangian density is the following: 

V{a,n) = imV- 0+^(0-0)' (3.62) 

for A fixed in the large N limit. Here is an 0{N) vector, = (cr, vf) and vr 
represents the — 1 pions. In what follows, we will consider two different cases of 
the potential y(cr, vf), with (m^ < 0) or without (m^ > 0) symmetry breaking. 

Our first order of business is to identify the correct order parameter for the 
phase transition and then to obtain its equation of motion. Let us define the 
fluctuation field operator x{^, t) as 

a = Mt)+xix,t), (3.63) 
with 00 (t) a c- number field defined by: 



. 1 f^.,M0mi, (3.64) 

nJ Trp(t) ^ ^ 

We can use the Liouville equation as before to arrive at the following evolution 
equation for the order parameter 0o(t): 
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To proceed further we have to determine the density matrix. Since the Liouville 
equation is first order in time we need only specify p(t = 0). At this stage we could 
proceed to a perturbative description of the dynamics (in a loop expansion). 

However, as we learned previously in a similar situation [|TB|, and as we have 
argued in this chapter, the non-equilibrium dynamics of the phase transition cannot 
be studied within perturbation theory. 

In the presence of a vacuum expectation value, the Hartree factorization is 
somewhat subtle. We will make a series of assumptions that we feel are quite 
reasonable but which, of course, may fail to hold under some circumstances and 
for which we do not have an a priori justification. These are the following: i) no 
cross correlations between the pions and the sigma field, and ii) that the two point 
correlation functions of the pions are diagonal in isospin space, where by isospin 
we now refer to the unbroken 0{N) symmetry under which the pions transform 
as a triplet. These assumptions lead to the following Hartree factorization of the 
non-linear terms in the Hamiltonian: 



x' - 


-> 6{x^)x^ + constant 


(3.66) 


x' - 




(3.67) 


7? ■ vf)^ - 


(^2 + — — ^ (vf^)??^ + constant 


(3.68) 


->2 2 


- {^')x' + ^'{x') 


(3.69) 






(3.70) 



where by "constant" we mean the operator independent expectation values of the 
composite operators which will not enter into the time evolution equation of the 
order parameter. 

In this approximation the resulting Hamiltonian is quadratic, with a linear term 
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in X- 



Hnit) 



+ 



Kit) 2 , Mlit) 



TT 



(3.71) 



Here 11^, 11^ are the canonical momenta conjugate to x{^)^ ^(^) respectively 
and is recognized as the derivative of the Hartree "effective potential" p5| , Q 
with respect to 00 (it is the derivative of the non-gradient terms of the effective 



action H, [T8|, |25[). 

To obtain a large N limit, we define 

Mt) = 

with 



(3.72) 
(3.73) 



{^') ^ 0{1) ■ ix') ^ 0{1) ; <P ^ 0{1). (3.74) 

We will approximate further by neglecting the 0{j^) terms in the formal large N 
limit. We now obtain 



Mlit) 
Mlit) 



N<p{t) 



A 



A 

2' 
A 



(3.75) 

(3.76) 
(3.77) 



We obtain the following set of equations 
0(t) + 0(t) 



(27r)3 2^1,(0) 



(3.78) 
(3.79) 



92 



The initial conditions for the mode functions U^{t) are 

f/+(0) = l ; U^{0) = -tuj^,{0) (3.80) 

Since in this approximation, the dynamics for the tt and x fields decouples, and 
the dynamics of x does not infiuence that of or the mode functions and {ip"^) we 
will only concentrate on the solution for the vf fields. 

3.4.1 Renormalization 

The renormalization procedure is exactly the same as that for the Hartree case in 
the previous section (see section 3.3.1). We carry out the same renormalization 
prescription and subtraction at t = as in the last section. Thus we find the 
following equations of motion 



2 



R 



t/+(t)= 0(3.81) 



with the initial conditions given by eq.( |3.80| ) and with the subtracted expectation 
value given in section 3.3.1. 

In contrast with the Hartree equations in the previous section, the cutoff de- 
pendence in the term proportional to (p^ in eq.( p.81| ) has disappeared. This is a 
consequence of the large N limit and the Ward identities, which are now obvious 
at the level of the renormalized equations of motion. There is still a very weak 
cutoff dependence because of the triviality issue which is not relieved in the large 
N limit, but again, this theory only makes sense as a low-energy cutoff theory. 
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3.4.2 Numerical Analysis 

Unbroken symmetry 



To solve the evolution equations numerically, we now introduce dimensionless quan- 



tities as in eq. (|3.55| ) obtaining the following dimensionless equations: 



Tj + Tj + rf + grjT,{T) = 
+ q^ + l + r]^T)+gE{T) 



U^{t) = 



f/+(0) = 1 ; /7+(0) = -z^q^ + l + vm 



dr 



(3.82) 

(3.83) 



S(r) 



A/M, 

q dq 



X 



Kir] 



1, M 
+ - In 

g2 + 1 + 7^2(0) 2 \Mn 



(3.84) 



For particle production in the 0{N) model, the final expression of the expec- 
tation value of the number operator for each pion field in terms of dimensionless 
quantities is the same as in eq.( p.59|) , but the mode function U^{t) obeys the 
differential equation in eq. (|3.56|) together with the self consistent condition. 

Figures (lO.a-c), (ll.a-c) and (12.a-c) show r/(r) ; S(r) ; N{t) for r/(0) = 1;2 
g = 0.1; 0.3 and A/M/j = 100 (although again we did not find cutoff sensitivity). 
The dynamics is very similar to that of the single scalar field in the Hartree approx- 
imation, which is not surprising, since the equations are very similar (save for the 
coefficient of the cubic term in the equation for the field expectation value). Thus 
the analysis presented previously for the Hartree approximation remains valid in 
this case. 
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Broken symmetry 

The broken symmetry case corresponds to choosing = — < 0. As in the 
case of the Hartree approximation, we now choose fin as the scale to define dimen- 
sionless variables and renormalization scale. The equations of motion for the field 
expectation value and the mode functions now become 



■^7] -r] + rf + gr]J:{T) = 



f/+(r)=0 (3.85) 



with S(r) given by eq. (|3.84|) . As in the Hartree case, there is a subtlety with the 



boundary conditions for the mode functions because the presence of the instabilities 
at r = for the band of wavevectors < < 1, for ^^^(0) < 1. Following the 
discussion in the Hartree case (broken symmetry) we chose the initial conditions 
for the mode functions as 

d 

dr 



U^{0) = 1; —f/,+ (0) = -Vg2 + 1+7^2(0) (3.86) 



These boundary conditions correspond to preparing a gaussian state centered at 
?7(0) at r = and letting this initial state evolve in time in the "broken symmetry 
potential" [gg. 

Figures (13.a-c) show the dynamics of ?7(r), S(r) and N{r) for //(O) = 0.5 
g = 0.1 and cutoff A/fifi = 100. Strong damping behavior is evident, and the time 
scale of damping is correlated with the time scale for growth of S(r) and iV(r). 
The asymptotic value of rj^r) and S(r), 77(00) and S(oo) respectively satisfy 

- 1 + r/2(oo) +^S(cx)) = (3.87) 

as we confirmed numerically. Thus the mode functions are "massless" describing 
Goldstone bosons. Notice that this value also corresponds to V'{(f)) = in eq.( p.75| ). 
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An equilibrium self-consistent solution of the equations of motion for the field 
expectation value and the fluctuations is reached for r = oo. 

Figure (13. d) shows the number of particles produced per correlation volume 
as a function of (dimensionless) wavevector q a,t r — 200. We see that it is strongly 
peaked at g = clearly showing that the particle production mechanism is most 
efficient for long-wavelength Goldstone bosons. At long times, the contributions 
from g 7^ becomes small. Figures (14.a-c) show the evolution of r]{T), E(r) and 
N[t) for ri{0) = 10~"^, g = 10^^. The initial value of r] is very close to the "top" 
of the potential. Due to the small coupling and the small initial value of r]{0), 
the unstable modes (those for which < 1) grow for a long time making the 
fluctuations very large. However the fluctuation term S(r) is multiplied by a very 
small coupling and it has to grow for a long time to overcome the instabilities. 
During this time the field expectation value rolls down the potential hill, following 
a trajectory very close to the classical one. The classical turning point of the 
trajectory beginning very near the top of the potential hill, is close to rjtp — y/2. 
Notice that tjIt) exhibits a turning point (maximum) at r) ~ 0.45. Thus the 
turning point of the effective evolution equations is much closer to the origin. This 
phenomenon shows that the effective (non-local) potential is shallower than the 
classical potential, with the minimum moving closer to the origin as a function of 
time. 

If the energy for the field expectation value was absolutely conserved, the expec- 
tation value of the scalar field would bounce back to the initial point and oscillate 
between the two classical turning points. However, because the fluctuations are 
growing and energy is transferred to them from the q — mode, rj is slowed down 
as it bounces back, and tends to settle at a value very close to the origin (asymp- 
totically about 0.015). Figure (14. b) shows that the fluctuations grow initially and 
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stabilize at a value for which 5'S(oo) ^ 1. The period of explosive growth of the 
fluctuations is correlated with the strong oscillations at the maximum of ?7(r). This 
is the time when the fluctuations begin to effectively absorb the energy transferred 
by the field expectation value and when the damping mechanism begins to work. 
Again the asymptotic solution is such that —1 + i]'^{oo) + g'E{oo) = 0, and the 
particles produced are indeed Goldstone bosons but the value of the scalar field 
in the broken symmetry minimum is very small (classically it would be rjmin — 1, 
yet dynamically, the field settles at a value ^^(oo) ~ 0.015!! for g = 10~^). Figure 
(14.c) shows copious particle production, and the asymptotic final state is a highly 
excited state with a large number (C(10^)) of Goldstone bosons per correlation 
volume. 

The conclusion that we reach from the numerical analysis is that Goldstone 
bosons are extremely effective for dissipation and damping. Most of the initial 

potential energy of the field is converted into particles (Goldstone bosons) and the 
field expectation value comes to rest at long times at a position very close to the 
origin. 

Notice that the difference with the situation depicted in figures (13.a-c) is in the 
initial conditions and the strength of the coupling. In the case of stronger coupling, 
the fluctuations grow only for a short time because 5'S(t) becomes 0{1) in short 
time, dissipation begins to act rather rapidly and the expectation value rolls down 
only for a short span and comes to rest at a minimum of the effective action, having 
transferred all of its potential energy difference to produce Goldstone bosons. 

Figures (15.a-c) show a very dramatic picture. In this case ri{0) = 10"'^, g = 
10~^^. Now the the fleld begins very close to the top of the potential hill. This 
initial condition corresponds to a "slow-roll" scenario. The fluctuations must grow 
for a long time before 5'S(r) becomes 0{1), during which the field expectation 
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value evolves classically reaching the classical turning point at rjtp — -\/2, and then 
bouncing back. But by the time it gets near the origin again, the fluctuations have 
grown dramatically absorbing most of the energy of the field expectation value and 
completely damping its motion. In this case almost all the initial potential energy 
has been converted into particles. This is a remarkable result. The conclusion 
of this analysis is that the strong dissipation by Goldstone bosons dramatically 
changes the dynamics of the phase transition. For slow-roll initial conditions the 
scalar fieki relaxes to a final value which is very close to the origin. This is the 
minimum of the effective action, rather than the minimum of the tree-level effective 
potential. Thus dissipative effects by Goldstone bosons introduce a very strong 
dynamical correction of the effective action leading to a very shallow effective 
potential (the effective action for constant field configuration). The condition for 
this situation to happen is that the period of the classical trajectory is of the same 
order of magnitude as the time scale of growth for the fluctuations. 

The weak coupling estimate for this dynamical non-equilibrium time scale is 
Ts ln(l/5r)/2, which is obtained by requiring that gT,{T) ^ 1. For weak couphng 
the mode functions grow as U^{t) ~ e'^ for < 1, and S(r) ~ e^^. The numerical 
analysis conflrms this time scale for weak coupling. For weakly coupled theories, 
this non-equilibrium time scale is much larger than the static correlation length 
(in units in which c = 1) and the only relevant time scale for the dynamics. 

Clearly the outcome of the non-equilibrium evolution will depend on the initial 
conditions of the fleld expectation value. 

Our results pose a fundamental question: how is it possible to reconcile damping 
and dissipative behavior, as found in this work with time reversal invariance? 

In fact we see no contradiction for the following reason: the dynamics is com- 
pletely determined by the set of equations of motion for the fleld expectation value 
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and the mode functions for the fluctuations described above. These equations 
are solved by providing non equilibrium initial conditions on the field expectation 
value, its derivative and the mode functions and their derivatives at the initial 
time t — 0. The problem is then evolved in time by solving the coupled second 
order differential equations. We emphasize the fact that the equations are second 
order in time because these are time reversal invariant. Now consider evolving 
this set of equations up to a positive time to, at which we stop the integration 
and find the value of the field expectation value, its derivative, the value of all 
the mode functions and their derivatives. Because this is a system of differential 
equations which is second order in time, we can take these values at to as initial 
conditions at this time and evolve backwards in time reaching the initial values at 
t = 0. Notice that doing this involves beginning at a time to in an excited state 
with (generally) a large number of particles. The conditions at this particular time 
are such that the energy stored in this excited state is focused in the back reaction 
to the field expectation value that acquires this energy and whose amplitude will 
begin to grow. 

It is at this point where one recognizes the fundamental necessity of the "in-in" 
formalism in which the equations of motion are real and causal. 

This concludes our discussion of the deatiled numerical analysis of the dissipa- 
tion in self-interacting scalar field theories. Let us now turn to the conclusion of 
this chapter where we will recapitulate what we have achieved, place our work in 
perspective and discuss the future directions that need to be explored. 
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3.5 Conclusion 



In this chapter we have used the real time formahsm to study one of the most im- 
portant time dependent and non-equihbrium phenomena that follows every phase 
transition namely dissipation and reheating. Phase transitions play an important 
role in cosmology. The universe we live in is an interesting place because of out 
of equilibrium phenomena. [§] Thus, it is important to systematically analyze the 
non-equilibrium phenomena associated with phase transitions. 

The real time formalism of finite temperaure field theory is the natural tool to 
use in understanding the non-equilibrium phenomena associated with phase tran- 
sitions. We have previously used this formalism to study domain formation and 
growth, The dynamics of dissipation and reheating that occur during a phase 
transition are among the most interesting non-equilibrium phenomena. In particu- 
lar, investigating this in depth is important for understanding the consequences of 
inflation]^ and other cosmological phase transitions. This chapter was devoted to 
understanding the damping of coherent field oscillations due to the phenomenon 
of particle production which is the microscopic mechanism for quantum dissipa- 
tion. We studied the simplest case of self-interacting scalar filed theories to enable 
us to concentrate on the basics. Here, we considered a system with a discrete as 
well as another system with a continuous internal symmetry in both the unbro- 
ken and broken symmetry case. Our results show that damping of oscillations 
is most efficient for the broken continuous internal symmetry case because of the 
availibility of the goldstone channel for dissipation. 



Our investigations of non-equilibrium phenomena in phase transitions! 16 1 have 



consistently shown the neccesity of accounting for non-perturbative effects when 



studying the long-time behaviour of the system. |]20|] If one restricts oneself to the 
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one-loop evolution equations one quickly sees that the first order effects grow to 
become of the same order as the zeroth order effects thus invalidating the use of 
perturbation theory for studying late time effects. To extract the non-perturbative 
results we have used the self-consistent Hartree approximation. 

The Hartree approximation is implemented by factorising all powers of the 
field higher than two and expressing them in terms of products of two point cor- 
relation functions and upto quadratic terms in the field and demanding that a 
self-consistency condition be satisfied. Once this is done one still needs to take 
care of divergences in the integrals that appear in the evolution equations. This is 
done by renormalisation using the WKB technique. 

We must point out that the Hartree approximation is uncontrolled in the single 
scalar theory. This approximation becomes exact in the N ^ oo limit of an 0{N) 
model which we have discussed. In our investigations, we also take the large N 
limit that is known as a systematic and controlled expansion in powers of for 
which the subleading corrections can be studied. 

By solving,plotting and studying the solutions to the evolution equations we 
were able to understand many interesting features of interacting systems. In par- 
ticular, we gained an insight into the microscopic mechanisms of dissipation for 
interacting systems. The solution to the coupled integro-differential equations de- 
scribing the time evolution for these interacting system was obtained numerically. 

We have discussed two self-consistent non-perturbative approximation schemes 
in this paper. The Hartree approximation and the large N approximation for the 
0(N) model. Of these the large N approximation has the advantage that it is a 
controlled approximation to which the corrections can be calculated in a systematic 
manner. In fact, the solutions to the evolution equations obtained in the large N 
approximation are in fact also simpler to analyse and understand. In particular, 
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the phase relationships between the various quantities such as the zero mode, the 
fluctuations and the particle production can be picked up at a glance from the 
plots diplaying the solutions. 

By examining these plots we noticed that the frequency of oscillations of both 
the fluctuations and the particle production is double the frequency of the zero 
mode oscillations. You may do this simply by counting the peaks of the oscillations 
in a fixed interval of time for the various quantities. Further, notice that the peak 
in particle production takes place exactly when the zero mode is at the minimum 
of its potential and therefore has its maximum kinetic energy within its oscillation. 
Thus, particle production and dissipation are proportional to the magnitude of the 
velocity within each oscillation of the zero mode. Within one complete oscillation 
of the zero mode it passes through the lowest point in its potential twice and there 
are two peaks in particle production within each complete oscillation of the zero 
mode. 

This clear phase relationship between the zero mode and the particle production 
seen in the large N approximation gets somewhat more complicated for the cutoff 
dependent Hartree case because higher harmonics of the fundamental frequency 
start playing a significant role in the time dependence. 

However, in both cases if one realises that there is a fundamental time scale in 
the problem related to the oscillation time scale of the zero mode, then one can 
gain some further insight into the behavior of the system and draw some further 
conclusions. 

Thus, if one chooses to call the time-scale of one zero mode oscillation as the 
microscopic time-scale of the problem, one can average over this microscopic time- 
scale to get the behavior on much larger time-scales - the macroscopic behavior 
of the system. By re-examining the plots of particle production in this light one 
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can see that on macroscopic time-scales the time- averaged particle production is 
a monotonically growing function of time which reaches an asymptotic final value 
when all the modes are sufficiently populated. 

All the momentum modes are not equally populated. In fact, energetically 
one expects lower momentum modes to be more populated than higher momen- 
tum modes. By examining the plots of N(k) versus k, one can draw the following 
conclusions. One notices that there are some distinct resonances due to favor- 
able phase relationships at specific momenta compared to neighboring momenta 
in phase space. Further, the lowest energy resonance is always the dominant res- 
onance because the fields are bosonic fields and energetically it is preferable to 
occupy the low momentum modes. 

Further, the phase space factor k^dk makes the contribution of the low mo- 
mentum modes to the total occupation number less significant than the sum of 
occupation numbers over the higher momentum modes. This can be verified by 
examining the magnitudes of the numbers for N[k) and comparing it to the total 
occupation number obtained by integrating over the momenta, Ntot- 

In all the different cases examined by us as time progresses, the amplitude of the 
zero mode oscillations decreases as the fluctuations grow and particle production 
increases.The amount of particle production and damping increase as the contri- 
bution of the cubic term in the zero mode evolution equation becomes significant 
and also as the coupling is increased. 

We examined the behavior of the system in both the unbroken symmetry case as 
well as the broken symmetry case. In the 0(N) model in the large N approximation 
one can see a very striking difference between the unbroken symmetry case and 
the broken symmetry case. Whereas in the unbroken symmetry case there is some 
diminishing of the amplitude of the zero mode oscillations due to dissipation, the 
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zero mode oscillations never completely die out. However, for the broken symmetry 
case there is a very rapid damping of the oscillations of the zero mode about its 
central value. The amplitude of oscillations about this central value goes almost 
completely to zero very quickly. 

Whereas for moderate couplings the zero mode settles asymptotically to a non- 
zero final value in the broken symmetry case there is a further interesting twist in 
the extremely weak coupling case. Thus for instance for g ~ 10"^ what happens 
is that fluctuations grow by such a large amount in the time that the zero mode 
is still above the spinodal point that after making a brief excursion to non-zero 
values the zero mode returns to a zero value. This happens because it is driven 
there by the fluctuations term in the evolution equation. 

Thus, we see in all cases that the zero mode damps out as particle production 
takes place. Further, in the broken symmetry case, the final value of the zero mode 
depends on the coupling and is different from the position at the bottom of the 
tree level potential. 

Here then is the summary of all the results gleaned from the plots of the 
solutions to the evolution equations. The fundamental time-scale of the problem is 
set by the oscillation time-scale of the zero mode when it is placed in the potential 
and allowed to evolve. Within each such oscillation, particle production takes 
place prefentially when the zero mode is at its lowest position in the potential and 
thus has the greatest magnitude for its velocity and kinetic energy. Further, if one 
keeps in mind this fundamental time-scale in the problem and factors out all the 
oscillations in various quantities on this time-scale one can follow with the eye the 
mean value of the quantities if one is interested in the coarse-grained behaviour of 
the system. Thus, in all cases examined, the fluctuations grow, particle production 
increases and the zero mode amplitude decrease as time progresses until they reach 
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their final values. The distribution of occupation number in momentum space 
shows the presence of resonances. Further, the lowest energy resonance is always 
the dominant resonance. 

Our study also reconciles dissipation in the time evolution for the coarse grained 
variable with time reversal invariance, as the evolution is completely specified by 
an infinite set of ordinary second order differential equations in time with proper 
boundary conditions. 

Our formalism and techniques are sufficiently powerful to give a great deal of 
insight into the particulars of the dissipation process. In particular, we can see 
that the damping of the field expectation value ends as the particle production 
ends. This shows that our interpretation of the damping as being due to particle 
production is accurate. 

It is useful to compare what we have done here with other work on this issue. 



We have already mentioned the work of Calzetta, IIu|3^ and Paz[0. These 
authors use the closed time path formalism to arrive at the effective equations of 
motion for the expectation value of the field. Then they solve the perturbative 
equations and find dissipative evolution at short times. In particular Paz finds the 
kernel that we have found for the effective equations of motion of the expectation 
value in the perturbative and Hartree case. However, his perturbative solution is 
not consistent. 

We remedy this situation by studing the non-perturbative Hartree equations, 
which must necesarily be solved numerically as we do. 

There has been other, previous, work on the reheating problem, most notably 



by Abbott, Farhi and Wise0, Ringwald||36|| and Morikawa and Sasaki[^]. In all of 



these works, the standard effective action is used, so that the expectation value is 
of the "in-out" type and hence the equations are non-causal and contain imaginary 
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parts. In essence, they "find" dissipational behavior by adding an imaginary part 
to the frequency that appears in the mode equations. We see from our work (as well 
as that of Calzetta, Hu and Paz) that this is not necessary; dissipation can occur 
even when the system is evolving unitarily. This comment deserves a definition of 
what we call dissipation here: it is the energy transfer from the expectation value 
of the g = mode of the scalar field to the quantum fluctuations (g 7^ 0) resulting 
in damped evolution for the q = mode. 

To what extent are we truly treating the reheating problem of inflationary 
models? As stated in the introduction, reheating typically entails the decay of the 
infiaton into lighter particles during its oscillations. What we do here is understand 
how the quantum fluctuations and the ensuing particle production influence the 
dynamics of the evolution of the expectation value of the field. Thus technically 
speaking, this is not the reheating problem. However, we are able to understand 
where dissipation comes from in a field theory, and are able to give a quantitative 
description of the damping process for the expectation value of the field. 

The techniques we develop here are easily adapted to the case where the infiaton 
couples to fermions and also to other scalars, a case that we have also studied p4|. 
In this connection, during the course of this work, two related pieces of work on the 
reheating problem have appeared^]. They both look at the effect of particle pro- 
duction from the oscillations of the infiaton field due to parametric amplification. 
What they do not do is to account for the back reaction of the produced particles 
on the evolution of the expectation value of the infiaton. As we have learned with 
our study, this back reaction will eventually shut off the particle production, so 
that these authors may have overestimated the amount of particle production. 

So far our in our analysis of dissipation we have been working in Minkowski 
space. We must extend our analysis to an expanding universe. This is perhaps the 
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most interestinng and promising future direction in the study of the dissipation 
of vacuum energy following a phase transition. Some of the machinery for this 
extension has already been developed|T^. We now need to do a detailed numerical 
analysis of the process of dissipation via particle production in curved spacetime. 

Through the detailed applications presented in this chapter and in earlier works 
we hope we have taken another step in making finite temperature field theory in 
both the imaginary time and real time formalisms a part of the toolkit of the 
working cosmologist. Applications of the real time formalism to the study of other 



time dependent phenomena such as thermal activation] 35 1 and domain formation 



and growth 1 16 1 are also available for the reader who wihses to see more detailed 
examples. In this connection, extension of these studies to curved space-time, 
in particular, for domain formation and growth is again a useful and interesting 
direction for future investigations. 
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Figure Captions 

Note: In all figures that involve time recall that r is the dimensionless time 
defined in this chapter as r = niRt 

Fig.l: Damping of the oscillations of the field about the minimum of its poten- 
tial 

Fig. 2: First order quantum correction for discrete symmetry case r^i(T) for 
771(0) = ; 7)1(0) = ; rici(0) = 1 ; 7)^(0) = 0. The cutoff is A/uir = 100. 

Fig. 3: Tree level potential in the unbroken symmetry case. (For purposes of 
illustration we have shown the potential: V{(f)) — ^ + ^ ) 

Fig. 4: Tree level potential in the broken symmetry case. (For purposes of 
illustration we have shown the potential: V{(j)) = + ^ ) 

Fig. 5. a: 77 (r) vs r in the Hartree approximation, unbroken symmetry case. 
^ = 0.1 ; 77(0) = 1 ; A/Mr = 100. 

Fig. 5.b: S(r) vs r for the same case as in fig. (5. a). 

Fig. 5.c: N{r) vs r for the same case as in fig. (5. a). 

Fig. 6. a: tjIt) vs t in the Hartree approximation, unbroken symmetry case. 
^ = 0.1; 7^(0) = 4; A/M^ = 100. 

Fig. 6.b: E(r) vs r for the same case as in fig. (6. a). 
Fig. 6.c: N(t) vs t for the same case as in fig. (6. a). 

Fig. 7. a: 77(r) vs r in the Hartree approximation, unbroken symmetry case. 
^ = 0.1 ; 77(0) - 5 ; A/Mr = 100. 

Fig. 7.b: E(t) vs t for the same case as in fig. (7. a). 
Fig. 7.c: N{t) vs r for the same case as in fig. (7.a). 

Fig. 8. a: r){T) vs r in the Hartree approximation, unbroken symmetry case. 
g = 0.05 ; 7^(0) = 5 ; A/Mr = 100. 

Fig. 8.b: E(t) vs r for the same case as in fig. (8. a). 
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Fig. 8.c: N{t) vs t for the same case as in fig. (8. a). 

Fig. 9. a: ?7(r) vs r in the Hartree approximation, broken symmetry case. 
g = 10-5 ; 7^(0) = 10-5 ; Af/iR = 100. 

Fig. 9.b: E(t) vs t for the same case as in fig. (9. a). 
Fig. 9.c: N[t) vs r for the same case as in fig. (9. a). 

Fig. 10. a: rj{r) vs r in the large N approximation in the 0(N) model, unbroken 
symmetry case, g ^0.1 ; r]{0) = 1 ; A/Mr = 100. 

Fig. 10. b: S(r) vs r for the same case as in fig. (10. a). 
Fig. lO.c: N{r) vs r for the same case as in fig. (10. a). 

Fig. 11. a: r]{T) vs r in the large N approximation in the 0(N) model, unbroken 
symmetry case, g = 0.1 ; rj^O) = 2 ; A/Mr = 100. 

Fig. ll.b: E(t) vs t for the same case as in fig. (11. a). 
Fig. ll.c: N{t) vs t for the same case as in fig. (11. a). 

Fig. 12. a: ?7(r) vs r in the large N approximation in the 0(N) model, unbroken 
symmetry case, g ^ 0.3 ; r){0) = 1 ; A/Mr = 100. 

Fig. 12.b: E(t) vs t for the same case as in fig. (12. a). 
Fig. 12. c: N{t) vs r for the same case as in fig. (12. a). 

Fig. 13. a: rj{r) vs r in the large N approximation in the 0(N) model, broken 
symmetry case, g — 0.1 ; r}{0) — 0.5 ; A/hr — 100. 

Fig. 13. b: S(r) vs r for the same case as in fig. (13. a). 
Fig. 13.c: N{t) vs t for the same case as in fig. (13. a). 

Fig. 13.d: Number of particles in (dimensionless) wavevector q, Nq{T) at r = 
200. 

Fig. 14. a: rj{r) vs r in the large N approximation in the 0(N) model, broken 
symmetry case, g — 10~^ ; r}{0) — 10"^ ; A/hr — 100. 
Fig. 14. b: S(r) vs r for the same case as in fig. (14. a). 
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Fig. 14.c: N{t) vs t for the same case as in fig. (14. a). 

Fig. 15. a: ?7(r) vs r in the large N approximation in the 0(N) model, broken 
symmetry case, g = 10~^^ ; 7^(0) = 10"^ ; Af/iR = 100. 
Fig. 15.b: E(t) vs r for the same case as in fig. (15. a). 
Fig. 15. c: N{t) vs r for the same case as in fig. (15. a). 
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